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Abstract: We study the problem of moduli stabilisation in explicit heterotic orbifold 
compactifications, whose spectra contain the MSSM plus some vector-like exotics that can 
be decoupled. Considering all the bulk moduli, we obtain the 4D low energy effective action 
for the compactification, which has contributions from various, computable, perturbative 
and non-perturbative effects. Hidden sector gaugino condensation and string worldsheet 
instantons result in a combination of racetrack, KKLT-like and cusp-form contributions to 
the superpotential, which lift all the bulk moduli directions. We point out the properties 
observed in our concrete models, which tend to be missed when only "generic" features of 
a model are assumed. We search for interesting vacua and find several de Sitter solutions, 
but so far, they all turn out to be unstable. 
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1. Introduction 

In searching for the Standard Model (SM) of particle physics and its supersymmetric gen- 
eralisations within String Theory, a myriad of constructions have been explored. Promising 
models have been built using heterotic orbifold compactifications, intersecting D-branes or 
D-branes at singularities and F-theory GUTs, to name a few. What all these models have 
in common, is that progress was made by focussing on the realisation of phenomenologically 
viable gauge groups and chiral matter spectra and, in particular, by neglecting gravity and 
the dynamics of the compactification. 

Eventually, however, it will be necessary to embed the supersymmetric standard model 
(SSM) into globally consistent and stable string theory constructions. For example, super- 
symmetry breaking is usually associated with the dynamics of the moduli, and therefore 
in order to make predictions for the soft masses we require an understanding of how the 
moduli are stabilised. Ideas abound also in how to address the moduli stabilisation prob- 
lem, and some progress has also been made towards dove-tailing these constructions with 
those of the SSM, particularly in Type IIB models [1]. So far though, the possibility of 
realising the SSM in a genuine metastable string compactification seems out of reach. 

In this paper we continue towards this objective within the context of heterotic orb- 
ifold compactifications. Heterotic orbifolds are arguably amongst the simplest frameworks 
for such studies. They have had a long history as candidates for realistic string theory 
constructions, particularly in the eighties and nineties, and as such, much effort was also 
put towards understanding supersymmetry breaking and moduli stabilisation there. Re- 
cent years have seen a significant revival in this class of models, since it was realised that 
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the idea of orbifold GUTs applied to heterotic orbifolds makes an effective tool to uncover 
the MSSM in the latter [2]. 

This progress is marked by the studies of a particularly fertile patch in the heterotic 
landscape, in which around 300 models were found containing the MSSM spectrum, to- 
gether with only some vector- like exotics that can be decoupled [3]. They are built us- 
ing a Z^-II orbifold compactification, with a non-standard gauge embedding (leading to 
N = (0, 2) worldsheet supersymmetry) and some number of discrete Wilson lines. We shall 
consider the problem of moduli stabilisation in some of these concrete promising models. 

In their totality, these models contain, as well as the MSSM and the supergravity 
multiplets, a hidden gauge sector, a number of vector-like exotics, SM singlets (typically 
charged at least under Abelian factors of the hidden gauge sector), and the bulk moduli. 
Amongst the twisted fields there are also flat directions, and some of these may have 
a geometric interpretation as blow-up modes. The latter are known as twisted moduli, 
but, interestingly for the possibility that our universe is an orbifold, they are not always 
present [4]. As a first step, our focus will be on all the bulk moduli. Therefore, although we 
will begin by treating all degrees of freedom, at the end we shall assume that the vector-like 
exotics and SM singlets are decoupled at some scale close to the string scale, consistently 
with supersymmetry, e.g. due to dynamics such as those discussed in [3]. 

In order to study the bulk moduli dynamics, it is first necessary to construct the low 
energy effective field theory for the given orbifold compactification. Fortunately, there is 
a good knowledge of the effective Lagrangian for orbifold compactifications (for some nice 
reviews see [5,6]). The relevant quantities can be computed using dimensional reduction 
and conformal field theory techniques [7, 8]. Another powerful tool [9-11] is modular 
invariance, which in the simplest models corresponds to an SL{2, Z) symmetry for each 
geometric modulus, and which is expected to hold even non-perturbatively. In what follows 
we shall compute the leading terms in the low energy effective field theory that describes 
concrete MSSM orbifold compactifications. Then, building on the work of the eighties and 
nineties, since all directions can be lifted, it seems natural to expect that (meta)stable 
vacua can exist. 

Indeed, although the dilaton and geometric moduli are all flat directions at tree level, 
non-perturbative effects tend to lift them all. Since the dilaton describes the tree level 
gauge couplings, gaugino condensation in a hidden non- Abelian sector leads to a non-trivial 
potential for that field [12]. A racetrack potential [13-15] or KKLT-type potential [16-18] 
may then be sufficient to stabilise it. The gauge couplings also receive threshold corrections 
from massive string states that depend, via certain cusp-forms, on several of the Kahler 
moduli and all the complex structure moduli [8,19], and thus a gaugino condensate may 
also stabilise those fields or even force compactification [15,20-22]. Yukawa couplings 
between twisted fields arise thanks to worldsheet instantons [7], and are suppressed with 
the Kahler moduli that describe their separation. This lifts all the remaining Kahler 
moduli directions [23] . The most complete and general analysis to date on some of these 
effects can be found in the seminal work [15], which showed the success of racetrack dilaton 
stabilisation in the presence of hidden matter, together with the stabilisation of a universal 
Kahler modulus. Another interesting explicit analysis was performed in [22], where it was 
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shown that a single gaugino condensate is sufficient to stabilise several geometrical moduli 
via a cusp-form mechanism, although the vacuum expectation value (vev) of the dilaton 
was set by hand. 

The present work is a first attempt at treating all the bulk moduli together in explicit 
models that give rise to the MSSM. Our models have three Kiihler moduli, one complex 
structure modulus plus the dilaton, giving a total of 10 real degrees of freedom. The various 
ingredients just described, give rise to a superpotential for all the fields, composed of a 
mixture of three types of terms, which have often been used in three popular stabilisation 
mechanisms: racetrack, KKLT and cusp-form stabilisation. 

Several differences occur when considering concrete models compared to the inspiring 
toy models of the past, however. Some of these are: (%) the terms in the action are almost 
completely determined with very few free parameters (which come mainly from the little 
studied physics that decouples exotics), (ii) the modular symmetry groups are generically 
broken from SL(2,Z) to some congruence subgroup due to the presence of Wilson lines 
[24,25], (Hi) it is typically difficult to find more than one condensing gauge group, especially 
without decoupling hidden matter (since the latter tends to destroy asymptotic freedom), 
and subsequently, it is hard (but not impossible) to find dilatonic racetrack models, (iv) the 
moduli-dependent threshold corrections to the gauge couplings often do not appear with 
the required sign to force compactification, (v) it is not justified to take a universal Kahler 
modulus, and moreover the dynamics of all the bulk moduli are highly coupled together. 
All these features make the search for a metastable (de Sitter) vacuum more challenging 
than previously thought and our results underline this. 

Due to the complexity of the system, our search for minima is largely numerical. We 
do find many de Sitter vacua, but all of them have unstable directions. We study their 
properties, including the presence of almost flat directions. Whether or not a metastable 
de Sitter vacuum might be waiting to be found is currently under investigation. 

The organisation of the paper is as follows. In the next section, we review briefly the 
orbifold construction in heterotic strings. Expert orbifolders can skip this introduction, 
and go directly to Section [|, where we describe the low energy effective theory arising from 
generic heterotic orbifolds, which takes the form of an = 1 supergravity. In particular, 
we recount the various perturbative and non-perturbative effects that can contribute to the 
dynamics of the moduli. Building on this discussion, Section || is devoted to the detailed 
study of the problem of moduli stabilisation in two concrete 2^-11 orbifold models, which 
contain the MSSM. We present explicit de Sitter vacua. Finally in Section [5] we discuss 
our results and the prospects to resolve the important open questions. 

Throughout the paper we take K4 = 1 and use dimensionless quantities. Whenever 
units are necessary, we say so explicitly. 

2. An Invitation to Heterotic Orbifolds 

In this section we review briefly the compactification of the EsxEs heterotic string on 
non-freely acting, symmetric Tjn toroidal orbifolds [5,6], which are the setting for our 
investigation. 
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An orbifold is obtained by modding out a discrete set of isometries P of a 6D torus 
T 6 , characterised by a lattice A. P is called the point group and we assume the torus T 6 
to be factorisable, i.e. T e = T 2 x T 2 x T 2 . In general, the basis vectors e a of the lattice 
A can be chosen to be the simple roots of a Lie group, whose metric is g a p = e a ■ e@. The 
action of the point group generator •& on A is given by 

<& : e a -> ep @p a , (2.1) 

where O is the Coxeter element of the corresponding Lie algebra. In terms of the lattice 
metric, Eq. ( |2.1| ) becomes 

i? : 9ap -> Q T g a pQ ■ (2.2) 

An admissible lattice should be invariant under P. This restricts the choice of A, but leaves 
some free parameters, which correspond to the geometric moduli of the compactification. 
In 'En orbifolds, the generator d can be made diagonal, such that 

= diag(e 2 ™ 1 , e 2 ™ 2 , e 2 ™ 3 ) 

where v = (0, v\, V2, V3) is the twist vector that parametrises the action of the orbifold on 
the compact space . Since 1? is of order N, it follows that N Vi = mod 1. Additionally, 
v is subject to the constraint ^ V{ = that guarantees the holonomy to be a (discrete) 
subgroup of SU(3) (but larger than SU(2)) and consequently, M = 1 in 4D. 

In the bosonic formulation, points of the torus are given by six real or, equivalently, 
three complex coordinates Zj, % = 1,2,3. Under 1?, the torus coordinates transform as 

Zl A ti&i . (2.3) 

In order to include the fact that Z is a point on T 6 , i.e. Z = Z + A with A € A, it is 
customary to define the space group S as the semidirect product S = P x A. An arbitrary 
space group element g = (i? fc , n a e a ) acts then on the coordinates of the compact space as 

Z i) k Z + n a e a , n a £Z. (2.4) 

It is easy to see that the action of a non-trivial space group element leaves invariant the 
points / satisfying 

(1 - $ k ) f = n a e a , n Q €Z. (2.5) 

Such points are called fixed points. Sometimes, instead of /, one uses the associated space 
group element g to refer to the fixed points. Since all elements from the conjugacy class of 
g, defined as {hgh^ 1 \ h € S}, have an equivalent action on the compact space, then the 
corresponding fixed points are equivalent in the orbifold. 

For each point group element "& k , one obtains a different number of inequivalent fixed 
points. Therefore, usually one refers collectively to the set of fixed points for a given k as 
the d k twisted sector. For k = 0, the twist action is trivial and therefore there are no fixed 

1 Here vq = represents the trivial action of the twist on the Minkowski 4D space with complex index 
/i = in light-cone gauge coordinates. 
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points. This sector is called the untwisted sector. The action of the point group on the $ k 
twisted sector is conveniently encoded in the local twist vector k = kv. In orbifolds with 
non-prime N, there are some sectors in which the twist acts trivially on one of the three 
two-torii. In these cases, we get fixed torii instead of fixed points. The holonomy in these 
sectors is a subgroup of SU(2), which leaves N = 2 supersymmetry unbroken. 

Modular invariance requires the action of the space group to be accompanied by a 
corresponding action on the 16 gauge degrees of freedom Xj of the heterotic string. The 
action of $ can be described by a shift 2 

X^X + nV, (2.6) 

where V is the shift vector. The lattice translations e a are translated to the gauge degrees 
of freedom as discrete Wilson lines A Q . Thus, an arbitrary space group element g = 
($ fc , n a e a ) affects the gauge degrees of freedom as 

X ^ X + 7t(kV + n a A a ). (2.7) 

Not every shift and Wilson lines are compatible with the space group. First, from = 1 
it follows that NV must be a vector in the weight lattice of EsxEs, i.e. the shift has to 
be of order N too. Also the Wilson lines are constrained. Since the space group elements 
g = (1, n a e a ) and "d l gi)~ l = (1, -d e n a e a ) = (1, mpep) are in the same conjugacy class, 
the Wilson lines n a A a and tjirAr must coincide. Direct consequences of this are that not 
all the Wilson lines are non-trivial nor independent, they can take only quantised values, 
and their order N a is restricted by the geometry of the compact space. 

Demanding modular invariance of the one-loop vacuum-to- vacuum amplitudes imposes 
additional constraints on the orbifold parameters which, for TL^ orbifolds without Wilson 
lines, read [26] 

N(V 2 -v 2 ) = mod 2. (2.8) 

In the presence of non-trivial Wilson lines, Eq. ( |2.8| ) has to be amended by [27] (no sum- 
mation implied) 

N a A a ■ V = mod 2 , (2.9a) 
N a A 2 a = mod 2 , (2.9b) 
gcd(iV a , Np) A a -A p = mod 2 («//?), (2.9c) 

where gcd stands for the greatest common denominator. These conditions guarantee 
anomaly freedom of the low energy effective theory. 

The simplest gauge embedding of the orbifold action corresponds to identifying the shift 
vector V with the twist vector v and setting A a = or, in the context of supergravity, 
to embed the spin connection into the gauge connection. This is the trivial solution to 
the modular invariance condition fl2.8|) and it is called standard embedding (SE). On the 
string theory side, the standard embedding preserves (2, 2) worldsheet supersymmetry. 

2 Another admissible embedding is a 16D rotation in the gauge degrees of freedom. 
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From the phenomenological point of view, the SE is not very attractive, since it is not 
possible to get realistic spectra 3 . Therefore, we work with non-trivial solutions to fl2.8j), 
or non-standard embedding (NSE) models. On the string theory side, these correspond to 
(0, 2) supersymmetric theories on the worldsheet. Moreover, also from phenomenological 
requirements, we are interested in models with non-trivial discrete Wilson lines. Indeed, all 
orbifold models with promising phenomenological properties correspond to this category [3]. 
We describe below the massless spectrum of this type of models. 

2.1 Massless orbifold spectrum 

The massless spectrum of an orbifold is comprised of two kinds of closed strings. The 
untwisted states stem from the original spectrum of closed strings of the heterotic string. 
In the supergravity limit, they correspond to the components of the 10D supergravity 
multiplet and the EsxEs vector multiplets that are left invariant by the orbifold action. 
The twisted states arise from closed strings attached to the fixed points. These strings exist 
only due to the orbifold, i.e. they cannot be identified as part of the supergravity limit of the 
heterotic string. It is natural to associate each state with a space group element commonly 
called constructing element. E.g. states in the untwisted sector are "constructed" by the 
element (1, 0) whereas those lying at the fixed point of the origin in the k = 2 twisted 
sector are associated with (i9 2 , 0). 

The general form of an orbifold state associated with g = n a e a ) can be written as 

W) = \q s h)R ® P \Psh)L ® |/) • (2.10) 

In this formalism, q s ^ = q + k encodes the quantum numbers of the (right-moving) string in 
the physical degrees of freedom of the compact and Minkowski space, where q is an SO (8) 
weight (e.g. for k = q is a weight of either 8^ or 8 S ) and k = kv is the local twist vector 
defined above. /3 denotes a product of (left-moving) oscillator excitations /3 I _ 1 in the 16 
gauge dimensions, in the Minkowski space, ^1. and /31 1+ „. in the holomorphic i and 
antiholomorphic i directions of the compact space, respectively. Here r/j = fcj mod 1, such 
that < T)i < 1. Further, p s ^ = p + kV + n a A a (with p a weight of Eg xEg) is a 16D vector 
that carries the information on the gauge quantum numbers. Finally, the component |/) 
contains the geometric information regarding the localisation in the compact space of the 
orbifold state. 

The quantum numbers of a massless state must fulfill 

Z m l = hl+N -l + 5c = 0, 

\™\ = ksh - \ + = o , 

where 8c = ^ Yl% r ]i{^-~ r 1i) corresponds to a change in the zero point energy due to (twisted) 
oscillators. The oscillator number N is given by 

3 

n = Y,^ N l + ^ N T ( 2 - 12 ) 

i=i 



3 There are two reasons for this: first, the resulting 4D gauge group is Ee, much larger than the SM 
gauge group, and secondly the number of 27's is always larger than three. 
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with f\i = —ki mod 1, such that < fji < 1. Here, N l g and N* 1 count the number of 
excitations related to the constructing element g in the i and i directions. 

Massless states \ip) with constructing element g = ($ fc , n a e a ) acquire in general a 
non-trivial phase under the action of an arbitrary space group element h = m a e a ), 

h ^ g27ri [p Bh -(lV+m a A a )-'R-lv~^ ((kV+n a A a )-(lV+rn a A a )-k-lv)] |^ ^ 

where we have defined the vector of IZ-charges 

K = q sh -N g + N* . (2.14) 

The orbifold projection consists of letting space group elements h such that [g, h] = 
act on the states \ip). Only invariant states under this action are retained in the orbifold 
spectrum. 

2.1.1 Untwisted sector 

In the untwisted sector, k = implies that the masslessness condition ( |2.11| ) reduces to 

1 9 ~ 1 9 1 , 

^p 2 +N-l = = -q 2 --. (2.15) 

Therefore, the solutions to this equation are given by q 2 = 1 and p 2 = 2, N = or 
p 2 = 0, N = 1. The resulting massless states can be gathered together in the following 
categories: 

Gravity multiplet and dilaton. The states of the form \q)R <S> /3^ 1 |0)l (/u runs in the 
Minkowski coordinates) such that q ■ v = mod 1 correspond to a) the degrees of freedom 
of the graviton g^ v , b) the dilaton </>, and c) one surviving component of the 10D B 2- form, 
whose dual is the so-called model independent axion, a. Since the orbifold preserves M = 1 
supersymmetry, also the superpartners are included in this notation. We can build an 
J\f = 1 chiral supermultiplet from the dilaton and axion, whose bosonic component is given 
by 

S = ^det g e~ 2 ^ + ia (2.16) 

with g the internal metric. 

Geometric moduli. The geometric moduli in the effective theory are the pure internal 
components of the 10D supergravity multiplet that survive the orbifold projection. In our 
notation, they are expressed as \q)n ®/9^ 1 |0)l and \q)R ®PLi\0)l with (q — N + N*) -v = 
mod 1. Traditionally, the states \q) have been denoted as \i) or \i) depending on the 
position and sign of the vectorial SO(8) weight q. For instance, |1) denotes q = (0, 1,0,0) 
and |2) denotes q = (0, 0, —1, 0). In this notation, states with \q) = \i) and antiholomorphic 
oscillator index j are the h ' Kahler moduli. Since we will be working in the case with one 
Kahler modulus per 2-torus, then i and j coincide and the moduli are denoted as TJ. States 
with \q) = \i) and holomorphic oscillator index j are the h 2,1 complex structure moduli U. 

A more useful notation to illustrate the geometric origin of these moduli is achieved 
through the torus metric. This can be rewritten as g a p = e a ■ ep = R a Rp cos tp Q p, with R a 
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being the "radius" of the cycle e a , and ip a p the angle between e a and ep. According to 
Eq. ( |2.2j ), imposing P invariance amounts to requiring 

® T g*pQ = g a( s ■ (2.17) 

This requirement fixes some of the parameters R a and (p a p . The remaining free parameters 
are combined to build the h 2,1 Kahler moduli and h ' complex structure moduli. To this 
purpose, let us define the metric of the i th T?~ plane as 

g (i) = ( ff2i-l,2i-l 92i-l,2i | ^ 2 
V 92i-l,2i g2i,2i J 

In terms of this metric, the complex structure and Kahler moduli of the i th plane are given, 
respectively, by 



Ui = + i 4| , T l = \ (Jd^gV) + W^J) . (2.19) 

9 u 9n x 



Since Ui contain only ratios of terms in the metric, they describe deformations in the shape 
of the 2-torii. Information about the overall size of the torii is carried by T,. 

4D gauge group. The states of the form \q)n ® \p) with q ■ v = mod 1 and p ■ V = 

= p ■ A a together with the 16 Cartan generators of EgxEg \q)ji <S> /3£ 1 |0)x, I = 1, . . . , 16, 
correspond to the gauge bosons of the unbroken 4D gauge group G^p. Notice that since 
PL\\0)l does not transform under the orbifold action, the rank of the gauge group is 
always preserved 4 . However, out of the 480 original charged gauge bosons (p ^ 0) only a 
subset of them survives the orbifold projection. Therefore, the gauge group is in general 
a subgroup of EgxEg. For instance, SE induces the breaking Eg — >• Eg x G^d, with Gid = 
SU(3), SU(2) x U(l) or U(l) 2 depending on the Zat orbifold. In non-standard embedding 
with Wilson lines, the breaking takes a more complex form. 

Chiral multiplets. Untwisted charged matter arises from the states \q)R® \p) where the 
orbifold projection is given by p ■ (IV + m a A a ) — q ■ Iv = mod 1 with arbitrary integers 

1 and m a . Their gauge transformation properties with respect to G^d are encoded in p. 

2.1.2 Twisted sector 

Zero modes of the d k twisted sectors (k ^ 0) are associated to the constructing elements 
g = (t? fc , n a e a ) and, therefore, attached to the corresponding fixed points. These are 
matter states that take the form fl2.10| ) and are constrained by the general masslessness 
condition ( |2.11 ). As before, their gauge transformation properties are encoded in the gauge 



momenta p s ^- Note that, since at the fixed points the gauge symmetry can be larger than 
the one in 4D, from a local perspective, at the fixed points, the orbifold-invariant states 
transform as complete multiplets of the local gauge symmetry. This observation has been 
useful to develop the concept of local GUTs [2,29,30]. 



4 The rank of the gauge group can nonetheless be reduced e.g. by turning on continuous Wilson lines [28]. 
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In the appendices we give the complete 4D massless spectrum for specific heterotic 
orbifold compactifications. For all models of interest, the spectrum is composed of the su- 
pergravity multiplet and the MSSM spectrum, together with some "hidden" gauge groups, 
hidden matter, a number of vector-like exotic matter multiplets and the bulk moduli. 

2.2 String selection rules 

Couplings among matter fields are constrained by the symmetries of the orbifold. All in- 
variant interactions with non-vanishing strength can be identified by using the so-called 
selection rules [7,31], which arise from those symmetries. In this discussion 5 , we con- 
sider the coupling c = f\ a A a , where the matter fields A a are characterised by the gauge 
momenta p s h, a , "^.-charges TZ a and constructing element g a . 

Gauge invariance. Since the gauge quantum numbers of all matter states are encoded 
in the gauge momenta p s h, c is gauge invariant as long as it fulfills 

J>sh,a = 0. (2.20) 

a 

7£-charge conservation. This selection rule of stringy origin arises from the discrete 
symmetries of the 6D compactification lattice A that are left unbroken under the orbifold 
action. In factorisable orbifolds, c preserves these symmetries if it satisfies 

J2'Jl i a = -l mod Ni , for all i , (2.21) 

a 

where iVj are the orders of the T 2 orbifold twists, i.e. Ni are the smallest positive integers 
such that NiVi = mod 1 (no summation over i). 

Space group selection rule. String interactions are only possible if the closed strings 
involved can join and form other closed strings. The constructing element g a not only 
characterises the fixed point f a of the matter state A a but also describes the boundary 
conditions of the associated strings. Based on this, 

H[g a ] = (t, 0) , (2.22) 

a 

where [g a ] denotes some element of the conjugacy class of g a . The rotational part of this 
equation reads Y\k a $ ka = 1 or analogously Yl a k a = mod N, which is called point group 
selection rule. 

For trilinear couplings, it is convenient to express the space group selection rule in 
terms of the fixed points f a 

(i - e ka )(o p f a - Ta ) + e ka (t - e fc ' 9 )(eV/3 - t p ) = (i - e ka+k ^Q r f^ , (2.23) 

where r a $ are arbitrary lattice vectors and the integers p,q,r lie between and N/2. 

In the appendices, we give the Yukawa couplings allowed by the above selection rules 
for specific models. 

5 The index a, used here to denote the various matter fields, should not be confused with the lattice 
vector indices, e a . 
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3. Low Energy Effective Supergravity Theory 



To study the dynamics of the moduli fields, we construct the low energy effective field 
theory, which describes the physics at energies well below the compactification scale. This 
is described by a 4D N = 1 supergravity theory and thus we must identify the corresponding 
Kahler potential, K, gauge kinetic functions, f a , and superpotential, W. 

For the purely untwisted fields of the orbifold compactification, the action can be 
identified simply by directly truncating the low energy 10D supergravity action describing 
the heterotic string. Things are a bit more complicated for twisted sectors, which emerge 
due to the compactification. The approach, pioneered in [7,8], is to construct an effective 
field theory that yields the same scattering amplitudes as the full string theory does in the 
low energy limit. Since orbifolds represent exactly solvable superconformal field theories, 
all quantities can, in principle, be computed. 

The low energy effective action receives corrections from two perturbative expansions: 
the string loop corrections, in e 2 ^, and the supergravity approximation, in a'. These 
corrections depend on the dilaton and volume moduli respectively. However, 4D TV = 
1 supersymmetry gives rise to powerful "non-renormalisation theorems". The tree-level 
superpotential cannot receive corrections from (string loop or a') perturbative effects, and 
is corrected only by non-perturbative effects [32]. Similarly, the gauge kinetic function 
receives string loop corrections only at one- loop, and no higher [33]. 

Finally, the low energy effective theory enjoys a further constraint from the discrete 
target-space modular invariance, which is inherited from the underlying conformal field 
theory. This provides a valuable check on all computations. 

In this section, we first describe in detail the modular symmetry, and then present 
the terms of the low energy effective action for the light fields in a generic orbifold com- 
pactification. We consider the untwisted moduli, S,Ti,U m , gauge sector and matter fields, 
A a . When relevant, we distinguish between charged matter fields (charged under the SM 
or hidden non-abelian groups), Qj C A a , and SM and non-abelian singlets (charged only 
under hidden U(l)s), $r C A a . 

3.1 Target space modular symmetry 

The spectrum of states in an orbifold string compactification is invariant under certain 
discrete transformations on the moduli, together with the winding numbers and momenta. 
In the simplest orbifold compactifications, without Wilson lines, the corresponding target 
space duality group is SL(2, Z) T (M) x [<S'L(2,Z)] [ / 2 ' 1) . Under this symmetry, the moduli, 
4>j = (Ti,U m ), and matter fields, A a , transform as follows [11,34] 

fa Trl a i' 6 i' c i' d i £ Z > a i d i ~ h i c i = 1 t 3 - 1 ) 

ILjCpj -f (lj 



A a ->■ MapAp lj (icjtffj + dj) 



(3.2) 



The modular weights = (n l a ,ff^) are given by 6 [11] 

for q sh ,i = 0,-1. 



n *=\ qs Y , ;."; '' s "-' ;;• ; cua.) 

and 



- 1 - q s h^ + K- N* 1 for g sh)i / 0, - 1 , 

for q sh , m = 0,-1 



q s h,m-N™ + N* g m for g sh , m ^ 0,-1, (3-3b) 

where the internal momenta q s h and oscillator numbers N g ,N* have been defined in Sec- 
tion 2.1. The matrices M Q( g are field-independent, and describe how the twisted fields with 
the same weights and charges transform amongst each other. Although in the supergravity 
basis this transformation is not unitary, it is so for the canonically normalised fields. 

Together with the sigma-model symmetries, the target space modular symmetry can 
be anomalous. Part of this anomaly can be cancelled by a Green-Schwarz mechanism, 
which implies that the dilaton S transforms at 1-loop, as [35] 

S ^ S ~^2 E 4? lo g( ic #; + dj), (3.4) 

3 

where 5 3 GS are the gauge group independent Green-Schwarz coefficients describing 1-loop 
mixing between the dilaton and the geometric moduli. The remainder of the anomaly is 
cancelled by threshold corrections due to massive string states, which we discuss below. 
The Kahler potential, K, transforms at both tree and 1-loop level as 

K^K+ l °g\ ic j</>j + dj\ 2 , (3.5) 

i 

and, in order to obtain invariant matter couplings, the superpotential has to transform (up 
to a field-independent phase) as [23] 

W 

W -> . (3.6) 

Modular invariance means that several terms in the action can be written in terms of 
modular forms, which transform covariantly with some weight k. The most common of 
these is the well-known Dedekind eta function 



ry(0) = e -^/ 12 Y[ (1 - e~ 2n ^) , (3.7) 

n=l 

which is a cusp-form of weight k = 1/2. This means that it vanishes at real infinity 7 and 
its modular transformation reads (up to a field-independent phase factor) 

rj{(f>) -> (ic^ + d) 1/2 ?70). (3.8) 



6 In Eqs. (3.E), we have assumed that the moduli Ti (U m ) he in the i th ijn th ) torus, which is the relevant 
case for us. 

7 We take the cusp for SL(2, Z) at real infinity. Mathematicians instead take it at imaginary infinity. 
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The presence of discrete Wilson line backgrounds breaks the modular symmetries down 
to the subgroup that leaves the Wilson line invariant [24,25]. Typically, the surviving mod- 
ular symmetries are congruence subgroups of SX(2,Z), which further restrict the integer 
parameters, a, b, c, d, of the transformations. They can be computed by solving the con- 
straints found in [25]. For instance, the groups that we encounter below are 

T (N) : c = mod N 
T°(N) : b = mod N 

Ti(N) :a,d = l mod N and c = mod N 

T^iV) : a,d = 1 mod N and b = mod N (3.9) 

with of course also a,b,c,d € % and ad — be = 1. It is easy to show 8 that, for instance, 
rj(N(j)) transforms covariantly under Tq(N) but not under the full modular group SL(2, Z). 
For r°(iV), the corresponding function is rj(cp/N). 

We close this subsection by noting that modular invariance represents a powerful means 
to check the low energy effective theory that one computes. Moreover, identifying mod- 
ular functions with the required weights and zero/singularity structure provides a way to 
infer the moduli-dependence of various terms in the action, when direct computations are 
lacking. 

In the remainder of this section, we introduce the different contributions to the low 
energy N = 1 4D effective theory below the compactification scale. When reviewing 
previous computations, we present only the final results and refer the reader to the original 
literature for further details. 

3.2 Moduli and matter kinetic terms 

The kinetic terms for the moduli and matter fields are encoded in the Kahler potential. 
For the twisted fields they have been computed to lowest order in the matter fields A a , 
with 

\A a \ 2 J] W>i + 0;) n "<l- (3-10) 

3 

Then, including one-loop effects, we have [11,15,22,23,42] 

3 

- log(^ + ^-)+ \A a \ 2 H (h + h)^. (3.11) 



K = - log 



Here n« are the modular weights as defined before. It is easy to show that ( |3.11 ) transforms 



in the required way (3.5) under the modular symmetry. The Kahler potential receives 
further perturbative corrections as well as non-perturbative ones. These, however, should 
be small compared with the leading contribution above. 



8 Hint: write iV^ as {J+j^ and use f 
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3.3 Matter interactions 

At tree-level, there may be couplings between bulk fields, bulk fields and twisted fields, or 
twisted fields located at the same fixed points. All these interactions can be understood 
within field theory, and are of order one. At the same time, recalling the stringy nature 
of the twisted fields, even though their center of masses are localised, they can stretch 
away from the corresponding fixed points. This leads to non-perturbative corrections to 
the interactions between twisted fields, which can be interpreted as world-sheet instantons 
and which are exponentially suppressed in the Kahler moduli that measure the area of the 
wrapped instantons. 

We now discuss in detail the trilinear Yukawa couplings, and their contributions to the 
superpotential. Our focus here is on twisted fields only, which have a non-trivial dependence 
of the Kahler moduli. Then we briefly comment on the aspects of higher order couplings 
that are relevant for our purposes. 

3.3.1 Yukawa couplings 

The computation of trilinear Yukawa couplings was started in [7,36] and completed in [37]. 
They were also explicitly computed for several orbifolds in [38]. 

Selection rules for permitted couplings were discussed in Section |2?2| . A generic coupling 
allowed between three twisted fields from the sectors , &~^ k+ ^), located respectively 

at fixed points (/„, fp, / 7 ), is given by 9 , 

Y aM = F *Pi9s VVaYI II ^V r fc.,i, ex P Hr I cotfc ^ + cotZ.TrrVil 2 ) , (3.12) 

where the local twists are k = kv and I = Iv, and J is the set of indices j labelling invariant 
planes, with kj, lj or kj + lj = (mod 1). The summation is over the components of u 
orthogonal to the planes j £ J and is due to contributions from fixed points equivalent 
to / a ,/3,7 thanks to lattice translations. Furthermore, F a p^ = ^Jl a lpl^ is a normalisation 
factor, which counts the degeneracy due to fixed points related to the f a by rotations, with 
l a being the number of fixed points within the fundamental torus domain that are related 
to f a by powers of 0. Also, V\ is the volume of the unit lattice orthogonal to the invariant 
planes, with V\ = ^/det g, and \u\ 2 =u t gu. Finally, 

Ui G [/^-/ Q -^ + r a + (l-G fe+i )(l-es cd ( fc -'))- 1 A] 4 , (3.13) 

r J T(K)T(U)Y{l-k^U) , U < K t + < 1 , 

iM \rv fe<)1 _i i; i<k i + h<2, 

where fcj — > k{ mod 1 so that < ki < 1, and r Q denotes vectors in A satisfying the space 
group constraint for non-vanishing trilinear couplings 10 

(i - e k )(f a - r«) + e fc (i - e l )(fp - Tp ) = (i - e fc+/ )/ 7 • (3.15) 



9 One can check that this expression is invariant under permutations of a, f3, 7, as the physics dictates. 
We have done so for our explicit couplings in the following section. 



Note that Eq. (3.15) differs from (2.23) by some factors Q p . The reason is that the rotational contri- 



butions to Y a p 1 have been absorbed in the coefficient F a p 1 . 
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Let us briefly comment on how this result meets with our intuition on the couplings. 
As one can see from the expression ( 3.12[ ), the coupling strengths depend on the orbifold, 



the fixed points and the precise coupling considered. Moreover, the dependence on the 
Tj moduli is also determined by the particulars of the coupling, which are encoded in the 
vector |u| 2 as well as in V\. One can also see from this expression that, when the three 
twisted fields are located at the same fixed point, the leading contribution is perturbative 
and of order one. Instead, when one or all three fields lie on different fixed points, the 
coupling is exponentially suppressed. We calculate these couplings for concrete examples 
in Section |||, and observe all these features directly. 

We are now ready to construct a holomorphic superpotential from the above Yukawa 
couplings. For this we note that we have set to zero possible backgrounds in the anti- 
symmetric tensor, £?2i-l,2ij which provide the holomorphic completion for the argument 
of the exponential in ( 3.1 2| ) . Recall also that Eq. ( 3.12 ) was computed for canonically 



normalised matter fields. The supergravity formula, relating these physical couplings to 
the superpotential is [39] (we suppress the k, I indices) 

Y ^ = " k/2 (K -kTk-W ' <316) 

where K is the tree-level Kahler potential for the (untwisted) moduli fields. Inverting this 
expression, we obtain the holomorphic superpotential 

W yuk = h a ^(Ti)A a A[sA y . (3.17) 



Finally we note that, according to the modular symmetry transformations for A a ( [3.2; ) 
and W (|3.6| ), the Yukawa couplings must transform as [23,34] 

K^{Ti) -> M'l M~l M-,\ h a ,^(Ti) Hiicjh + dj yi-<-^-^ ( (3.i 8 ) 

j 

that is, the couplings transform, up to a non-trivial weight factor, amongst themselves for 
fields with the same modular weights and charges. 



3.3.2 Higher order couplings 

Higher order couplings are in principle also computable [40]. Any allowed couplings of 
order N, will have the generic form in the superpotential (in this section we put back 
explicitly the units) 

Wi n t = M P h (Ti)A ai A a2 . . . A aN . (3.19) 

where h aia2 .^ aN (Ti) has the required modular transformation properties and the fields A ai 
are dimensionless in Planck units. 
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The higher order terms N > 3, are non-renormalisable. However, they may give rise 
to effective mass terms if the coupling features N — 2 non-Abelian singlet matter fields, 
which acquire a vev at some high scale. In this way, we may assume that all SM exotics 
and charged hidden matter, Qi, acquire an effective mass with the mass matrix taking the 
form 

M u = m p h ri ...r r u(Ti) <$n • • • *iy> • (3-20) 

r 

This may be important to allow for gaugino condensation, discussed below. Assuming 
singlet vevs that lead to a universal mass scale for the SM exotics and charged hidden 
matter, the decoupling mass can be defined as 

M d = detM 1/dimM , (3.21) 

and Md < M s < Mp (M s = a' -1 / 2 is the string scale) for singlet vevs slightly below M s . 
We discuss this in a little more detail at the end of this section. 



3.4 Gauge sector dynamics 

The final sector that we need is the gauge sector. We focus on the hidden gauge group, 
G = Y\ a G a - The corresponding gauge couplings are field-dependent, and encoded in the 
gauge kinetic functions. In many cases, particularly when any charged hidden matter is 
decoupled at some intermediate high scale, < M s , the gauge coupling becomes strong at 
some scale, A < Below A, the gauginos condense, leading to another (field theoretical) 
non-perturbative contribution to the superpotential. 

3.4.1 Gauge kinetic function 

The tree-level gauge kinetic function, ft ree = k a S, with k a the level of the Kac-Moody 
algebra, receives corrections at one-loop, due to massive charged particles running in the 
loop: 

f a =k a S + A^(T i ) + A^(T i ,U m ). (3.22) 

Here, A^ d (Tj) are the field theoretical threshold corrections due to the massive charged 
hidden matter, and A^f s (Ti,U m ) are stringy threshold corrections due to massive and 
winding string states. The field theory contributions take the form [23] 

A^T,) = -^log^, (3.23) 

where the beta function coefficient, 

b a = -3C(G a ) + ^2T(Q a j), (3.24) 
i 

includes the charged matter, whilst 6^ = —3C(G a ) is that of the pure Yang Mills theory, 
valid below M<j. Here, C(G a ) and T{Q a I ) are the quadratic Casimirs in the adjoint and 
matter representations, respectively. 



- 16 - 



The stringy thresholds turn out to receive contributions only from N = 2 sectors in 
the string spectrum. In the absence of Wilson lines, they were computed in 11 [19,42], and 
can be written in terms of the Dedekind eta function as follows 

A a (T h U n ) = -^S ^) M=2 log V (^f . (3.25) 

Here j runs over all the complex structure and Kahler moduli associated to N = 2 planes. 
The Pj refers to the subgroup of the point group P, which leaves unrotated the j th complex 
plane of T 6 . Then T e /Pj is an orbifold with M = 2 4D supersymmetry and (^)-^ =2 are 
the corresponding N = 2 beta function coefficients, given by 

(V a f= 2 = -2C(G a ) + £ T(R a a ) , (3.26) 

a 

where denote the representations w.r.t. the non-abelian group G a , and the summation 
runs over all half-hypermutiplets of the theory. We discuss the modifications of these 
formulae due to the presence of Wilson lines in the following section. 

Notice that, since the Dedekind function vanishes at infinity, the threshold corrections 
diverge there. This divergence is due to the fact the Kaluza-Klein modes all become 
massless as the torus volume goes to infinity, signalling a breakdown of the 4D effective 
theory. The modular symmetry ensures there is a similar divergence as the torus volume 
goes to zero, when the winding modes become massless. If there are special points in the 
moduli space where extra charged modes become massless, we expect additional, milder 
singularities in the thresholds at those points. 

Recalling that the stringy threshold corrections help to cancel the modular and a- model 
anomalies, we can relate their coefficients to those of the anomalies, b£ , in the following 
way [19,42] 



\3 

The anomaly coefficients are computed to be 



il (Kf= 2 = b'j-k a 5 GS . (3.27) 



b'j = -C(G a ) +£r(Qj)(l + 2n\) (3.28) 



and so the relation ( J3.27 ) provides a way to calculate 5 GS . 

Having established the holomorphic gauge kinetic function, the expression for the 
running, loop corrected gauge coupling constant g a at a scale fj, is [23] 

7* + log -# + log 



9l 167T2 6 V ^ ) 16VT 2 & \Ml 



16vr 2 



11 There is generically an additional universal contribution to the thresholds, which also depends on the 
moduli [19], but has been little studied (see [41]). We neglect it in the present work, but it would be 
interesting to consider its effects in further studies. 
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Y,( b 'a -kaS' GS ) log\ V (^)\\ (3.29) 



16vr 2 

where the gauge coupling at the string scale is given by 

J> = S + 5 " ^2 E d hs log (fc + fc) , (3-30) 

and again, the sum over j runs over h 2,1 moduli and Kahler moduli associated to J\f = 2 
planes. Notice that, whereas in the supergravity holomorphic gauge kinetic function the 
threshold corrections kick in at the fundamental Planck scale, for the real gauge couplings 
g a the relevant scale is the string scale. Moreover, the gauge couplings must not only 
be real, but also modular invariant. Modular invariance is guaranteed by the effect of 



massless gauginos running in the loop, which give rise to the fourth term in (3.29). For 
an enlightening discussion on holomorphic gauge kinetic functions and modular invariant 
gauge couplings, see [43]. 

3.4.2 Gaugino condensation 

In the effective theory below the decoupling scale, Md, we are left with pure Yang-Mills 
groups. In this case, the gauge couplings generically become strong at some scale A a and 
the following superpotential is generated 

W gc ttj^da exp (^ffa) , (3-31) 

where d a is a constant that arises in the process of integrating out the condensate, and is 
given, up to an unknown constant c, by 

c b° 

d a = (3.32) 
e96vr 2 v ; 

in Planck units (see e.g. [5]) and e is just the Euler number. 

It is a simple exercise to show how W gc transforms under modular transformations 
(|3.6|) . Using the above expressions for /„ fl3.22D -( j3~28|) , the definition of Md ( |3.21 ) and the 
modular transformation of det M that follows from fl3.2j ) and ( |3.18 ), 



detA^^detTW ]J ( ic j^ + d j)~ El(1+2n ^ , (3.33) 
j 

one can verify that W gc transforms in the required way ([O]), provided that the charged 
hidden matter is such that 

J2nQi)^[(l + 2n j J )}-dimM^2[(l + 2n j J )T(Q J )} =0. (3.34) 
I J J 

This appears to be a stringy constraint on the matter content, which explicit models - 
consistent with all our assumptions - should fulfill. In the models analysed in Section 
we show that this non-trivial condition is indeed satisfied. 
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3.5 The scalar potential 

With the Kahler potential, K, superpotential, W, and gauge kinetic functions, f a , in hand, 
we can now compute the scalar potential of the action. There are two kinds of contributions: 

V = V F + V D . (3.35) 

The F-term potential is given by 

V F = e K (k aB D a WD b W - ?>\W\ 2 \ , (3.36) 



where the indices A,B label all the chiral supermultiplets present and D^W = 8aW + 
OaKW. For our generic Kahler potential, Eq. ( |3.11 ), and a general W, Q3.36 ) takes the 
form 

2 



V F = e 



K 



\YW S -W\ 2 + Y, 



1 



W Aa +A a ]l(<i> i + $ i )<W 



+E 



Y 



Y — 5 i — Y Y^ a l^l 2 < Uj(4>d + <Pi) ni 

S { W s + {fa + &)W^ -W-J2 <A a W Aa 



3\W\' 



, (3.37) 



+ 



and 5 3 



where we have defined Y = S + S — J^j & log {^j ~r rjy ■ - -< 

Additionally, there are D-term contributions to the scalar potential from each gauge 
multiplet, abelian and non-abelian alike. In general, these can be written as 



V, 



Da 



KatfA^ 



(3.38) 



where ta^ are the generators of the gauge group, appropriate to the representation in 
which the chiral superfields A a lie. However, amongst all the gauge multiplets, there is 
one special one, which is associated to a pseudo-anomalous U(l)x gauge symmetry. The 
would-be anomaly is cancelled by a Green-Schwarz mechanism, which requires the dilaton 
S to transform non-linearly under the U(l)x at one-loop, and which is associated to a field- 
dependent Fayet-Iliopoulos term. In this case, the contribution to the D-term potential 
is 

2 



x 



K a q a A a + 



(3.39) 



Here, q a are the U(l)x charges, and Sqs the U(l)x Green-Schwarz anomaly coefficient, 
which is given by [44] 



1 



96vr 2 



trtx , 



(3.40) 



where tx is the generator of U(l)x- An important remark is in order here. It is clear 
from the expression above that all geometric moduli appear in the D-term potential via 
K a , unless the modular weights vanish. In orbifold models, matter fields always have non- 
vanishing modular weights in at least one torus, thus this dependence is always present. 
Moreover, a non-trivial moduli dependency is also ensured via the perturbative correction 
to the dilaton appearing in Y. 
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3.6 Decoupling of exotics 

To complete this section, we make a few more remarks on the decoupling of exotics, first 
introduced in Subsection 3.3.2. Recall that the orbifold compactification initially gives rise 
to a spectrum which includes not only the particles of the MSSM, hidden gauge group 
and moduli, but also a number of vector-like exotic matter and SM singlets, which may be 
charged under non-Abelian hidden sectors or only under hidden U(l)'s. These must all be 
decoupled somehow if we are to describe Nature. 

One way that this may happen is if non-trivial vevs are induced for some non-Abelian 
singlets, thanks to e.g. the Fayet-Iliopoulos term described above. As explained in Sub- 
section 3.3.2 , the singlet vevs subsequently give masses to the exotics and charged hidden 
matter, due to their higher order couplings. At the same time, the hidden U(l) gauge 
symmetries will typically all be broken. This mechanism has been studied in [3], where it 
was argued that the vevs can be turned on in a way that preserves super symmetry, along 
the lines of [45]. Important to our picture, is that the vevs induced by the D-terms can be 
arranged to be constant with respect to the moduli and that the coupling strengths that 
give rise to effective mass terms are also moduli independent. Indeed, otherwise, we would 
not have effective mass terms but couplings between the exotics and the light scalar moduli 
fields. Although such moduli independence is not generic, it does seem to be possible, and, 
moreover, an interesting way to allow for an additional scale in the problem. It would be 
important to investigate these issues further 12 . 

For the time being, we note some significant phenomenological motivations for our 
introduction of an additional scale. On one hand we would like the exotic matter to be 
decoupled at least at the GUT scale, in order to maintain the gauge coupling unification of 
the MSSM. On the other hand, we would like the supersymmetry breaking scale, usually 
associated with the moduli stabilisation, to be down at the electroweak scale in order 
to understand the latter 13 . Furthermore, the decoupling of charged hidden matter at 
some high scale can be important to ensure that gauginos condense later, since massless 
charged matter would drive the jV=l beta- function coefficients towards positive values 
which destroy asymptotic freedom. This especially applies to gauge groups of small rank, 
as they tend to come with a lot of hidden matter in heterotic orbifold compactifications. 
Indeed it would otherwise be very difficult to find multiple condensing gauge groups in 
the minilandscape. Finally, it should be remembered that gaugino condensation in the 
presence of massless charged hidden matter is actually poorly understood (see e.g. [15]). 

Therefore, in the following, we simply assume that all the exotics are decoupled below 
some scale near the string scale, Mj, consistently with supersymmetry. In our low energy 

12 The possibility to dynamically decouple matter fields via FI D-terms in our explicit models turns out 
to be highly involved. Note that, as mentioned before, the FI D-terms have a non-trivial dependence on 
all the moduli, via their modular weights and 1-loop corrections. This possibility certainly deserves further 
investigation. 

13 One alternative (but less motivated) picture is that the moduli are somehow stabilised at a high scale, 
the exotics are then decoupled via the singlet vevs somewhere below this scale as in [3], and supersymmetry 
is broken later by some so far unknown mechanism. As yet another possibility, albeit technically unfeasible, 
all the fields may be fixed at the same time. FI D-terms [46] and matter fields active during moduli 
stabilisation [46,47] have been noted to help produce de Sitter/Minkowski vacua. 
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4D effective supergravity theory the physics of the decoupling is then parameterised by the 
singlet vevs and the scale Indeed, since supersymmetry is preserved, we can simply 
integrate out the heavy matter by replacing them with their vevs in the Kahler potential, 
superpotential and gauge kinetic functions. The action that results then describes the 
MSSM sector and the bulk moduli, and (since the "hidden" U(l)'s are all broken) the 
relevant contributions to the scalar potential can only be from F-terms. 

4. Moduli Stabilisation in Realistic Z 6 II Orbifold Models 

We are now ready to study the dynamics of the moduli that results from the action de- 
scribed above. In this section, we compute explicitly all the supergravity ingredients that 
we reviewed in previous sections, in concrete orbifold models. We stress here that whilst 
most of the ingredients have been discussed in the past, this is the first time that such re- 
sults are applied to "realistic" orbifold examples in their full glory. Perhaps unsurprisingly, 
we find that several of the features that have been used in the toy-models constructed to 
date, tend not to be observed in real models. At the same time, we will see that whilst the 
scalar potential seems to have sufficient structure to stabilise all of the bulk moduli, we 
found only unstable vacua. Whether or not these instabilities are inevitable is currently 
under investigation [48]. Meanwhile, since all the fields typically contribute to the breaking 
of supersymmetry, we generically obtain de Sitter vacua. 

We start by making some general observations and then proceed by analysing explicit 
models in which the exact MSSM can be found. 

4.1 Moduli stabilising contributions 

Let us first explain how the ingredients discussed in the previous section give rise to a 
scalar potential for all the untwisted moduli. 

As described above, the Yukawa couplings between three twisted fields receive world- 
sheet non-perturbative instanton contributions, which depend on the area of the cycles 
that the instantons wrap. This gives rise to a potential for the Kahler moduli of the 
corresponding cycles. Notice that if any of the twisted sectors involved has an invariant 
plane, then the corresponding twisted fields are not localised in that plane, so there is 
no area suppression in the couplings there. Therefore, W yu k depends only on the Kahler 
moduli belonging to planes that are rotated by the twists involved. 

Meanwhile, it is well-known that gaugino condensates in hidden gauge groups give 
rise to a non-trivial potential for the dilaton, S. Moreover, the threshold corrections to 
the gauge couplings lead to a dependence in W gc on all the moduli describing 2-torii that 
remain invariant under one of the orbifold twists. In this way, the gaugino condensates 
provide a non-trivial potential for all the complex structure moduli and some of the Kahler 
moduli, in particular all those that did not appear in the Yukawa couplings. Therefore, 
in principle, none of the S,Ti,U m are actually flat directions. We now ask, in concrete 
models, if it is possible to stabilise all moduli and to obtain a de Sitter vacuum. 
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Figure 1: Geometry of the T 6 torus of a "Zq-II orbifold. The big dots denote points left 
invariant by the points denoted by a star and a square are fixed only under $ 2 and $ 3 
respectively. In addition, the shorter /blue arrows denote the non-trivial Wilson lines given 
in Appendices A and B. 



4.2 Explicit Z 6 II MSSM models 

In this section we consider two representative MSSM models from the Zq-II miniland- 
scape [3], which exhibit however different generic properties for the moduli dynamics. They 
serve as good examples of the sort of structure one expects to see in the fertile heterotic 
orbifolds 14 . We now introduce the essential features of the models with the details left for 
the appendices. 

The twist for the Zq-11 orbifold is v = |(0, 1, 2, —3), for which an admissible choice 
of A is the root lattice of G 2 x SU(3) x SO(4) depicted in Figure |]. The action of the 
corresponding Coxeter element on the lattice basis is given by 

6ei = 2ei + 3e 2 , 0e 2 = —e\ - e 2 , 
0e 3 = e 4 , 9e 4 = -e 3 - e 4 , 

©e 5 = -e 5 , 9e 6 = -e 6 . (4.1) 



Point group invariance of the metric, Eq. ( 2.17| ), implies that there are five free parame- 



ters (associated with three real Kahler moduli and one complex structure modulus). In total 
we have the following relations among the parameters (in string units, i.e. Rf = Rf/a') 

a/3 1 
R 1 = y/3R 2 , R 3 = R±, cosip 12 = — — , cos99 34 = -- (4.2) 

and the free parameters are R\, i? 3 , i?5, Rq, cos ip^. Subsequently, using the standard 
definition of Kahler and complex structure moduli, Eq. ( 2.19| ), we find that the geometric 
parts of our moduli are given by 

Re^) =t 1 = ^, Re(T 2 ) = t 2 = & R* , Re(T 3 ) = t 3 = ±i? 5 Re sin^e , 

U 3 = u 3 + i^ 3 = ^| sin <p 56 + i ^ cos ip 56 , (4.3) 
where ip$Q is the angle between e$ and e%. 



14 Although the models presented here are based on Eg xEg heterotic Z6-II orbifolds, similar constructions 
exist in the context of its SO(32) sister and other TLn orbifolds [49]. 
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The different models are obtained by choosing different gauge embeddings (i.e. different 
shift vectors and Wilson lines). We now present each in turn, and write down explicitly 
their low energy effective theories by applying the results of Section 3. We assume for 
simplicity that non-Abelian singlets acquire vevs as discussed in Subsection [T(], which are 
constant A a = (A a ) and of about the same order. Moreover, we assume that all the exotics 
and charged hidden matter are subsequently decoupled at the same scale, We then 
focus on the moduli and hidden gauge group only, since the observable matter must have 
vanishing expectation values in order to preserve the SM gauge group. Our dynamical 
fields are thus S,Ti,Us, i = 1,2,3. 

4.2.1 Model I: double gaugino condensate 

The shift vector and Wilson lines for this model are given in Appendix |A|. The unbroken 
gauge group after compactification is SU(3) x SU(2) x U(l)y x [SO(8) x SU(3) x U(l) 6 ]. 
The full 4D massless spectrum is given in Appendix |A.2| . As explained above, the presence 
of discrete Wilson lines breaks the modular symmetry down to a congruence subgroup of 
SL(2, Z). In the present model, we compute the modular symmetry group to be SL(2, Z) x 
Ti(3) x Tq (4) x r°(4), where the transformations act on Ti,T%,T3 and U respectively. We 



give some useful information about these groups in Appendix A. 4. The resulting action is 
as follows. 



Kahler potential. The Kahler potential is given by ( 3.11| ), where the modular weights 



can be found in Table |3| and the Green-Schwarz coefficients are computed in Appendix A. 7. 
They are 

*GS = -T' *GS = -h 6 GS = 0. (4.4) 

In principle, all the SM singlet vevs contribute to K. Since they must be small in order for 
us to maintain perturbative control, we can safely neglect them at this level. 

Gauge kinetic function. We can compute the gauge kinetic function by referring to 



Section 3.4.1. The values for the beta- function coefficients that measure the field theoretic 
threshold corrections turn out to be 15 

bso(8) = -9, i>so(8) = ~ 18 

W) = 1 > b° sm = -9. (4.5) 

Meanwhile, although stringy threshold corrections in the presence of discrete Wilson lines 
have been scarcely studied 16 , we can infer their modular dependence from the modular 
symmetries and expected singularity structure. For example, assuming that the thresholds 
depend on T2 as log[r/(3T2) 2 ] leads to the required covariance under the modular symmetry 
Ti(3) in the second plane, and reproduces the expected divergences as the volume of that 
plane goes to zero or infinity (see the paragraphs below Eq. ( |3.26| )) 17 . Moreover, as in the 



15 Here we used C(SO(N)) = N — 2, T(N so(iv) ) = 1, C(SU(N)) = N, T(N SU(N) ) = ±. 
16 See the second reference of [24], and very recently [50]. 

17 We do not exclude other dependencies, for instance one that is covariant under Fi(3) but not To(3). 
Indeed, direct computations of the stringy thresholds in the presence of Wilson lines would be very valuable. 
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case without Wilson lines, we expect the contributing sectors to be the N = 2 ones, and 
this allows us to compute the coefficients of the thresholds. We study the corresponding 



auxiliary N = 2 theories in Appendix A. 6 



Putting all information together, the final result for the two hidden gauge groups, 
SO (8) and SU(S), are respectively 

h = S - ^ log + A log 77 (3T 2 ) 2 - JL (log rj (4T 3 ) 2 + log V (U/A) 2 

5 i M d 1 , ,„^2 3 / , . , TTIA . 2 



/3 = S " 4^2 log M^" 12^2 log??(3r2) 8^ (l0g7?(4T3) + lo S^/ 4 ) ) • ( 46 ) 



At this point one can check explicitly, using the data in the Appendix |A.2j , that the 
constraint (3.34) is indeed satisfied, and the expressions are modular covariant. Another 



nice check of the threshold coefficients is to note that computations for all the gauge groups 
lead to the same Green-Schwarz anomaly coefficients for the sigma-model and modular 
anomalies, via Eq. (|3.27|) . 

Superpotential. The superpotential for the moduli receives contributions from Yukawa 
couplings, higher order matter couplings and gaugino condensation. In the moduli poten- 
tial, we neglect the higher order couplings, since they will be suppressed with respect to 
the trilinear Yukawa couplings by the Planck scale and the small matter vevs (see how- 
ever [51]). Also, we consider only couplings between (non-Abelian) singlets, since we allow 
non-trivial vevs only for those fields. 

The allowed Yukawa couplings are computed in Appendix |A.3| . All those between 
gauge singlets are of the kind $ 2 # 5 'i9 5 . From ( [3.12 )-( 2.23| ), we can explicitly compute this 



coupling in terms of the Tj moduli (see Appendix^ for details). The result is: 

Y a 2 | 7 = F aM g s VhhN 255 ^2 ex P T i " 6 «2 T 2) » ( 4 - 7 ) 



where 



2 14 / 3 7T 4 

N255 = 3W(l/3) ' (48) 
b Ul = ttu\ L\ u\ , b U2 = ttu 2 L 2 u 2 , (4.9) 





L '=(- 3 2 J' i2= (-l 2 J' <410) 

The 2D lattice vectors in depend on the fixed points (f a , fp, / 7 ) of the fields involved in 



the coupling, and are given in C.2, as are the normalisation factors F a ^. 

The corresponding contributions to the holomorphic superpotential are computed using 
( |4.7| ) and ( 3.16Q . In Appendix we explain the necessary tricks to show that the explicit 



expressions are modular covariant, as they must be. Here, we consider only the leading 
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order terms in the sum over instantons, and replace the SM singlets by their vevs A a . The 
result is 18 

W yuk = 2N 255 Aj e - 2 -^/3 ^ 2 + _ (4n) 

Finally, with the gauge kinetic functions in hand, we can immediately write down 
the contribution to the superpotential that arises when the gauginos in the hidden sector 
condense: 

M d \ 3/2 



--e3^ e ~ 87r2S/3 {wS" V(3T2)4/9 [T]m) ^ Um2 ' (412) 
where we have used the definition of d a in ( 3.32| ). Notice here that several of the Dedekind 
eta functions, originating from the stringy threshold corrections, appear with positive pow- 
ers, in contrast to what is usually assumed in the literature. We comment below on the 
implications of this for moduli stabilisation. 

Before looking at the moduli stabilisation and de Sitter vacua, we present a second 
representative model, which has a single hidden sector gauge group. 

4.2.2 Model II: single gaugino condensate 

The shift vector, Wilson lines, massless spectrum and all other details of this model 
are given in Appendix [b|. The unbroken gauge group after compactification is SU(3) x 
SU(2) x U(l)y x [SO(14) x U(l) 5 ] . The modular symmetry is broken by the Wilson lines 
to SL{2,'L) x Ti(3) x T\{2) x T 1 (2), where the transformations act on Ti,T 2 ,T% and U 
respectively. 

The computation of the action follows exactly as for Model I, and again we relegate 
the details to the appendix. 

Kahler potential. The Kahler potential is given by ( 3.11| ), where the modular weights 
can be found in Table ||, and the Green-Schwarz coefficients are 

Shs = -U, *GS = -y, ^5 = 0- (4.13) 
Again, we neglect contributions to K from the singlet vevs, since they are small. 

Gauge kinetic function. The beta-function coefficients that measure the field theoretic 
threshold corrections are 

&SO(u)= -29, b° so(u) = -36. (4.14) 

Thus, the gauge kinetic function for the single hidden sector is given by 

7 i M d , 7 , , nm . 2 , 5 / ,„_. 2 , , ,__.„, 2 



fu = S - ^ log j£ + ^ log r, (3 T 2 y + g-j (log r, (2 T 3 ) z + log V (U/2) z ) . (4.15) 

Using the data in Appendix |b|, we can check that the constraint ( 3.34| ) is satisfied, the 
expressions are modular covariant, and the universal 5q S are obtained. 

18 We have taken the prototype couplings ngns2n64, 71337152 rae4, n^n^n^, 712371597166 (see Table and 
assumed consistently, for simplicity, the vevs 7152 = ne,4 — 7159 = ne,a = Ai, Tig — A2, \/2n23 = ^3- 
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Superpotential. As in the previous model, we consider the contributions to the super- 
potential from Yukawa couplings and gaugino condensation. The allowed Yukawa couplings 
for this model are computed in Appendix |B|. As before, couplings among singlets are all 
between iT 2 -!? 5 ?? 5 sectors, and the expressions for the couplings in terms of the Tj moduli 
can be found in Appendix ^|. Using the notation in the appendix, we can have couplings 
of four types. 

We can compute the holomorphic superpotential and check that the expressions are 
modular covariant. In what follows, we consider only the following terms 19 



yuk 



2N- 



255 ^-1 



A? (A 2 + A 3 e- 2 ^ 3 



(4.16) 



where we have considered only the leading instanton contributions, and replaced the singlets 
by their vevs. Notice the important fact that a constant piece has naturally arisen, coming 
from the leading contribution to the Yukawa coupling between twisted fields localised at 



the same fixed point (see Appendix C.2 ). 

Finally, the gaugino condensation contribution to the superpotential is simply, 



c 3 



-2w 2 S/3 



j^) r?(3T 2 )- 14 / 9 [r,(2T 3 )r,(U/2)\ 



-5/6 



(4.17) 



Notice that in this case, the Dedekind eta functions do appear with negative powers in the 
superpotential, as is usually assumed. 



4.3 Towards stabilisation of moduli and de Sitter vacua 

In this section, we study the problem of moduli stabilisation and de Sitter vacua in the 
two explicit orbifold models discussed in the previous section. 

The scalar potential that governs the dynamics of the moduli, after having integrated 
out the hidden matter, is given to a good approximation (in Planck units) by 



V F 



K 



Y 



Ws 
W 



1 



+E 



Y 



Y -&> 



W ^3 Y3 > w 



1 



(4.18) 



T u T 2 ,T 3 ,U,6i 



ft. 



, and now 



where, as before, Y = S+S—^j & log {4>j + <j)j) with <f>j = ±x 
W = W yu k + W gc - We have to analyse this potential for the five complex fields S, Ti, T 2 , T3 
and U. Due to the complexity of the system, we can only make limited progress analytically, 
in particular we can solve for the axionic parts. Thereafter, our search for vacua will be 
numerical. First, however, we can make some general observations. 



4.3.1 de Sitter vacua 

Before considering the moduli stabilisation mechanisms that may be at work, let us begin 
by discussing the structure of the potential, and in particular the possibility of obtaining de 



19 These come from the prototype couplings riTnesneg, nnnesn^g, n7nsin62, nnnein^ (see Table 0), by 
taking consistently the vevs nei = = ne,$ — uqq = Ai, nr = Aq., \f2nv; = As. We have also turned on 
other allowed couplings, but the results were not more interesting. 
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Sitter solutions. From ( 4.18[) we see that the positivity of the potential depends basically 



on how large the ratios W a /W can be, compared to one. If three of our five fields have 
W a /W sufficiently smaller or larger than one, then the potential will certainly be positive. 
In order for it to be negative, the ratio W a /W has to be of order one for at least three of 
the fields. One can then analyse the dependence of the ratio W a /W for each field. For this, 
let us recall the orders of magnitude that we are considering, i) The vevs of the matter 
fields must be around A Q < 1CT 1 in order for the calculation of the Kahler potential to 
be valid, ii) The decoupling scale must be lower than the Planck scale and thus we have 
M d < 10- 1 M P . 

For Model II, it is simple to estimate that the contribution from the dilaton will 
generically give rise to a constant bigger than one, and thus can cancel a big part of the 
negative piece in the potential. Moreover, for T%, the ratio W a /W can be seen to be much 
less than one, thus contributing again to the positivity of V . For the other fields, it depends 
on rj and thus can in principle be of order one and therefore give a null contribution to the 
potential. The upshot is that for Model II we can generically expect V > 0, and thus, de 
Sitter vacua. An analogous analysis can be made for Model I, with similar results. Note 
that the loop contribution from 5 3 is always smaller than the others, as it must be. 

We now ask whether or not we can expect to find vacua with all the moduli stabilised. 

4.3.2 The modular form mechanism for T2,T^,U 

Let us first focus on the three moduli, T2,Tz,U, which appear in the superpotential due 



to stringy threshold corrections, via the Dedekind eta functions defined in (3.7). It was 
suggested in [20] that such a dependence may be enough to force compactification and 
moduli stabilisation. In particular, if the eta functions appear in the superpotential with 
negative powers, then the scalar potential diverges as the relevant moduli go to zero and in- 
finity (this can be seen by recalling the asymptotic behaviour of the eta function described 
in Section |3.ip . Then a minimum in these directions is guaranteed. Moreover, with the 
simplest factorisable superpotentials such as those corresponding to a single gaugino con- 
densate only, the fixed points under the modular transformations are necessarily extrema 
of the potential [20,22]. 

Although the latter property is not observed for our models, which have more compli- 
cated superpotentials, Model II does have the required asymptotic behaviour to ensure the 
existence of minima in the T2,T^,U directions. However, in Model I, which is also a very 
common one in the landscape, the Dedekind eta functions appear in the superpotential 
with positive powers, so that the scalar potential goes to zero as the relevant moduli go 
to zero and infinity. In that case, we cannot ensure the existence of minima. Notice also 
that in realistic orbifolds like the ones we consider here, Wilson lines break the modular 
symmetries in each plane in different ways, and thus generically, it is not correct to assume 
equality of all the Kahler moduli (the dependence of the superpotential on each Kahler 
modulus assumes a different functional form). 

The existence of multiple gaugino condensates and how the threshold corrections enter 
the superpotential are somehow correlated with each other. If the Eg x Eg gauge group 
is broken down to yield multiple non-Abelian hidden gauge groups, the latter will tend to 
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be small, and thus one can expect to find large amounts of hidden matter. At the same 
time, much hidden matter, whether massless or massive, drives the auxiliary M = 2 beta- 
function coefficients towards positive values, which in turn means that the eta functions 
appear with positive powers in W gc . 

The case of T 2 in Model I is slightly different since it also always appear in the Yukawa 
couplings. Thus a different mechanism may help to stabilise this field. 

4.3.3 Racetrack and KKLT for S,T 1 ,T 2 

The dependence of the potential on T\ can occur only via the non-perturbative Yukawa 
couplings in Z^-H orbifolds. The form of the superpotential is then reminiscent of the 
KKLT one [16], albeit with coefficients that depend on the other moduli. The same is true 
for S, in the case of a single gaugino condensate. Otherwise, S may be stabilised with 
a racetrack induced by multiple gaugino condensates (see [15]). Similarly, if T 2 appears 
in W yu k, it may have a racetrack behaviour, with W yu k competing against the leading 
dependence in W gc . 

Since the standard racetrack and KKLT stabilisation mechanisms both require some 
degree of tuning of the coefficients, and since in our case those coefficients are actually field- 
dependent, whether or not they are successful in producing a minimum in all directions 
must be studied explicitly. We now turn to that analysis. 

4.3.4 Stabilisation of the axions 

We can make significant progress by first considering the imaginary parts of our five complex 
fields, the axions. These turn out to have periodic potentials, with multiple minima and 
maxima, several of which can be found analytically. To begin with, recall that the axions 
contribute to the scalar potential only via the superpotential. 

For the axions entering the potential via the Dedekind eta functions, it is in general 
difficult to identify the positions of all the extrema. For example, only for the simplest 
superpotentials do the fixed points correspond to extrema. However, it turns out that 
Im.(j)j = are typically minima, and so we may choose those points. This simplifies the 
problem for the remaining fields 20 . For example, consider the superpotential for Model II, 
which can now be written in the form 

W = d e~ ClTl + C 2 (T 2 , T 3 , U) e- c * s + C 3 (4.19) 

with all the coefficients real. This is reminiscent of a KKLT superpotential [16] for the 
fields T\ and S, although here the KKLT coefficients are no longer constants, but depend 
on the other fields. For the dependence on T 2 ,T3 and U, let us consider just the leading 
terms in the eta functions. Then the superpotential is simply a sum of exponentials, and 
the derivatives that appear in the scalar potential are easy to compute. 

Similar to the basic KKLT model, the resulting scalar potential depends on Im T\ = a\ 
and ImS* = c 5 as terms proportional to cosciOi, cosc2(Js and cos(ci<7i — c 2 eg ) . The 

2(, These extrema are independent of the real components of tfij. It is however clear that there are other 
critical points Im <f>j which do depend on the specific values of Re <j>j ■ 
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magnitudes and signs of the coefficients depend on C\, C2 , C3 as well as the real parts 
of the five complex fields. In any case, it is clear that o\ = niir/ci and as = ri2ir/c2, 
with ni,ri2 integers, are critical points, and so they represent interesting candidates for 
the minima. Notice that analogous arguments can be made for models that fall into the 
racetrack class. 

With extrema for the five imaginary parts in hand, we now turn to a numerical analysis 
for the real parts. We have also scanned numerically for all ten real and imaginary parts 
together. 

4.3.5 Numerical results and explicit dS vacua 

We have tried several methods to search numerically for minima in the scalar potential 21 . 
The most successful one we found was to search for solutions to the critical point equa- 
tions 22 dV/dx A = for x A the real (and imaginary) parts of S,Ti, U. 

By varying the input parameters c, A a and we have identified several critical points 
consistent with our approximations 23 in both models. In Tables p] and [2] we show some 
representative vacua of Model I and Model II, respectively. Mathematica codes with all 
the details of these vacua can be found at our webpage [52]. 

The vacua we find are all de Sitter, however they all turn out to have at least one 
unstable direction in the moduli space. This is evident from the plot of the scalar potential 
as a function of ReT2 in Model I, and/or from the mass eigenvalues, displayed in Tables [l] 
and One might be tempted to think that a shift to the left in V(ReT2) of Model 
I leads to the true minimum. However, this is not sufficient since any change in ReT2 
affects the whole scalar potential, resulting in a vacuum with (more) instabilities in other 
directions. In general, we have found it very difficult to avoid instabilities in the direction 
Re (T2, T3, U). It would be important to understand the reason for these instabilities, and 
whether or not they are inevitable. 

An interesting feature of both models is that there are some almost flat directions, as 
observed from the hierarchically small masses. The corresponding eigenvectors are linear 
combinations of Im S,ImT^ and ImU (and also ImT2 for Model II). The reason for these 
almost flat directions is as follows. First, notice that in Model I, the contribution to W gc 
from the SU(3) gaugino condensate is typically highly suppressed with respect to the SO(8) 
condensate due to the stronger exponential suppression in S, and the stronger suppression 
in Md/Mp. Thus, effectively, we have only a single condensate in both models. Moreover, 
the higher order contributions to the eta functions in both Model I and Model II, at their 
solutions, are strongly suppressed for T3 and U (the same can be said for T2, when ReT2 

21 We especially thank Yashar Akrami, David Grellscheid, Stuart Raby and Timm Wrase for detailed 
discussions on this point. 

22 We have performed numerous checks to our solutions e.g. minimizing numerically the one dimensional 
problems, examining the 2D and 3D plots, and increasing the numerical precision to up to 1000 figures. 
Note that not all candidate critical points that one finds survive these checks. For example, it can happen 
that the initial numerical computation finds minima in two directions that are slightly shifted from each 
other, and does not resolve this difference. 

23 In addition, we have taken the first 20 terms in the expansion of the eta function. We have verified 
that higher order terms do not alter our results. 
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c= 1/10 


Parameters 
20 Ay = 100 A 2 = A 3 = 1/10 


M d = 1/65 


4>i 


Moduli vevs 
Re(^) 


lm(0i) 


Ti 


3.51166 


3/2 


T 2 


0.24201 


-1/3 


T 3 


7.05981 


35/8 


U 


112.95695 


-506 


S 


0.09385 


1079/144tt 


Cosmological constant A = 6.45 x 10 19 


Effective dilaton Y = 0.32262 


Pi ~ <t>i 


Mass eigenstates 




Pi 


1.87 x 10~ 18 


3.57 x 10- 18 


P2 


-4.76 x 10~ 17 


1.2 x 10" 17 


P3 


3.15 x 10~ 20 


3.29 x 10~ 33 


PU 


9.84 x 10~ 23 


1.51 x 10~ 94 


PS 


1.14 x 10" 16 


2.16 x 10" 16 


Tachyon: Rep 2 ~ 0.9ReT 2 + 0.4ReS 



Table 1: An unstable de Sitter solution of Model I. The dominant contribution to the mass 
eigenstate pi arises from the modulus 0j. In particular, the tachyon is dominated by the 
Re T2 direction, as can also be appreciated in Figure H We give the solution to 5 significant 
figures, but we have computed it to a precision of 1000. The value of A is given in Planck 
units. 



becomes larger than one as in Model II). Therefore, to a good approximation, ImS, ImT3 
and ImU (and ImT2 for Model II) appear in the superpotential, and thus in the scalar 
potential, only via a single linear combination, which, e.g. for Model I, can be written as 

47T^ 7T 7T 

X = ImS--ImT 3 ImU. (4.20) 

3 9 144 v ' 

As a consequence, to a good approximation, only this linear combination of the three 
(or four) fields is lifted by the non-perturbative dynamics, leaving two (or three) linearly 
independent combinations as flat directions, and two (or three) corresponding shift sym- 
metries. Of course, the higher order corrections do lift the latter directions, rendering 
them only almost flat. We expect this feature to be also present in any metastable minima 
that may exist, being quite a generic consequence of the moduli stabilisation via gaugino 
condensation with threshold corrections. Such almost flat directions could potentially be 
interesting, since extremely light scalar fields are often called upon in cosmological models, 
for inflation, quintessence and the like. 

Finally, we comment on the vevs of the moduli that we find. Of course, due to the 
instabilities, our vacua cannot be considered as candidates for our universe, and so we have 
not attempted to obtain realistic values for the moduli at this stage. However, consistency 
of our analysis does impose some constraints. In particular, the overall volume must be 
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V(S) 

Re(S) 
0.0938558 



6.45 




144 n 

Im(S) 

Figure 2: Scalar potental as a function of the bulk moduli of Model I for the vacuum 
presented in Table |]. The scalar potential is rescaled by a factor 10 19 . 



large enough for the supergravity approximation to hold, and the dilaton large enough for 
the string loop expansion to be a good one. The latter is difficult to achieve, as has long been 
known in heterotic orbifolds. Meanwhile, it is interesting that anistropic compactifications 
emerge quite naturally in our setup, since these have been proposed to explain the mild 
hierarchy between the GUT and string scales. 
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Parameters 






c= 1/2 


A 1 ,A 2 ,A 3 = 1/2 


M d = 1/600 


c= 1/200 


A X ,A 2 ,A 3 = 1/2 


M d = 10" 3 






Moduli 


vevs 






4>i 


Re(^) 


lm{cf>i) 


(pi 


Re(^) 


lm((j)i) 


Ti 


0.60565 


3/2 


Ti 


0.60213 


3/2 


T 2 


2.05917 


19/6 


T 2 


3.99948 


19/6 


T 3 


6.48601 


37/4 


T 3 


11.95957 


37/4 


U 


25.94402 


-3 


U 


47.83828 


-3 


S 


0.53034 


809/144tt 


S 


0.97758 


809/144tt 


Cosmological constant A 


= 7.41 x 10~ 6 


A = 5.12 x 10~ 7 


Effective dilaton Y 


= 1.15312 


Y = 2.07762 






Mass eigenstates 






Pi ~ 4>% 


m Re« 




Pi ~ 4>i 




2 

m lm Pi 


Pi 


2.88 x 10" 5 


2.28 x 10~ 4 


pi 


2.57 x 10~ 6 


2.02 x 10^ 5 


Pi 


-1.67 x 10~ 5 


8.34 x 10~ 19 


pi 


-1.92 x 10~ 7 


8.61 x 10~ 36 


P3 


5.57 x 10~ 7 


1.21 x 10~ 37 


p-i 


1.2 x 10~ 8 


1.23 x 10" 68 


Pu 


2.84 x 10~ 8 


7.57 x 10~ 39 


Pu 


5.97 x 10" 10 


7.72 x 10~ 70 


PS 


4.27 x 10~ 4 


7.87 x 10~ 4 


PS 


5.28 x 10~ 5 


6.72 x 10~ 5 


Tachyon 


Rep 2 ~ 0.9ReT 2 + 


0.3ReT 3 + OAReS 


Rep 2 ~ 


0.9ReT 2 + 0.4ReT 3 


+ O^RcS 1 



Table 2: Unstable de Sitter solutions to Model II. As in Model I, the dominant contribution 
to the tachyon arises from ReT2. The potential V(ReT2) differs from the one for Model 
I because the tachyon has some components in some other real directions. A change in 
the input parameters provides a larger value for Y. We give the solution to 5 significant 
figures, but we have computed it to a precision of 1000. 



5. Discussion 



We have revisited the problem of moduli stabilisation in heterotic orbifold compactifica- 
tions, in the light of the recent discovery of fertile patches in the heterotic landscape where 
the MSSM spectrum can be found. After reviewing the derivation of the low energy effec- 
tive action describing general orbifold compactifications, we computed it explicitly for two 
concrete MSSM models. Our resulting actions survived several non-trivial checks, such 
as modular invariance, including a new stringy constraint on the matter content in the 
presence of gaugino condensation ( 3.34 ). Finally, we studied the dynamics for all the bulk 
moduli fields, making use of numerical techniques. We found several de Sitter vacua, but 
unfortunately they all suffer from instabilities. Curiously, we did not uncover any anti-de 
Sitter vacua that were consistent with our assumptions. Since the scalar potential grows 
rapidly, we expect them to appear at very small moduli vevs, which are not consistent 
with our approximations. On the other hand, the values of the cosmological constant at 
the vacua tend to be very large, although we should recall that this is counting only the 
classical value for the vacuum energy. 

Given that we have found it difficult to locate metastable vacua, it is interesting to 
compare our results with previous toy models for heterotic orbifold compactifications. The 
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most recent of these is [18], which treated two bulk moduli, S and T, together with a 
complex matter field. Although they included a D-term contribution, it turns out to be 
vanishingly small, so that their setup is very similar to ours. Indeed, it is easy to find 
metastable minima to their F-term potential. Thus it seems that the challenge we meet 
comes from considering all five bulk moduli fields that are present and the comparatively 
small number of free parameters. 

Our analysis leads us to question if new ingredients, beyond those already considered 
in the literature, may well be necessary in order to stabilise all the moduli. Along those 
lines, recall that we have neglected some moduli-dependent contributions to the effective 
Lagrangian, which are known to be present, but have been little studied. Firstly, there is 
a universal contribution to the threshold corrections [41], and secondly an overall factor in 
the Yukawa couplings, which can depend on bulk moduli of the untwisted planes [37]. It 
would also be very useful to have in hand direct and explicit computations of the threshold 
corrections in the presence of discrete Wilson lines [50]. Other effects that depend on the 
moduli and merit further study are those arising from the Casimir Energy and Coleman- 
Weinberg potential, which become relevant after supersymmetry breaking (see e.g. [53]). 

Otherwise, the main unknown in our analysis has been the dynamics that decouples 
exotics and hidden matter, which we parametrised with some (non-Abelian) singlet vevs 
and the decoupling mass scale. The introduction of this high scale, which separates the 
dynamics of the matter and moduli, is important for several reasons. Phenomenologically, 
we know that the exotics must be decoupled at least at the GUT scale to maintain gauge 
coupling unification, and yet the gauge hierarchy problem suggests supersymmetry breaking 
should occur down at the TeV scale. Practically, it renders the system tractable; indeed 
already with 10 real bulk moduli the search for minima has proven to be a significant 
challenge. Occasionally, it is also important because hidden matter, if not massive, tends 
to destroy asymptotic freedom in the hidden gauge groups, preventing the gauginos from 
condensing there. 

We remark that in considering only the bulk moduli stabilisation, we have assumed 
that all the other flat directions, corresponding to twisted fields, may be lifted in the 
decoupling process, and this certainly deserves more attention. One can argue that, since 
the orbifold MSSM candidates arise from M = (0, 2) models, the massless twisted fields 
might not generically correspond to blow-up modes [4]; thus it should be possible to stabilise 
them and remain within the orbifold description where we have much more control. A 
detailed treatment of the decoupling dynamics would require, among other things, explicit 
computations of the higher order couplings between non-Abelian singlet fields, and non- 
Abelian singlets and charged matter. For progress in this direction see [40]. Notice that if 
a decoupling scale is in play, D-term potentials are no longer relevant for the bulk moduli 
dynamics, since all the U(l) gauge groups are broken when the non-Abelian singlets acquire 
vevs (all the non-Abelian singlets are generically charged under all the hidden U(l)'s). 

Although unstable, the dS vacua we have found do have some interesting properties. 
In particular, the non-perturbative dynamics (i.e. gaugino condensation with threshold 
corrections) tends to give rise to some almost flat directions, and consequently hierarchically 
light fields. We expect that this feature can also hold for any metastable dS vacua that 
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might exist. Light fields could be useful for cosmology, and it would then be important to 
understand how quantum corrections contribute to the light masses after supersymmetry 
breaking. 

So far we must say that our results on the (non-)existence of metastable de Sitter 
in MSSM heterotic orbifolds are inconclusive. We are now developing more sophisticated 
numerical codes to scan our scalar potentials for metastable minima, and are widening our 
search in the minilandscape of MSSM models [48]. Then we must also ask if a model can be 
found with realistic values for the moduli vevs, including Re S ~ 2, an appropriate gravitino 
mass and cosmological constant. If and when such a vacuum is found, it would have a host 
of applications. We could begin to study in detail the phenomenological and cosmological 
aspects of the compactification, starting with supersymmetry breaking and the low energy 
sparticle spectrum, models of inflation, resolutions to the cosmological moduli problem and 
so on. 
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A. Details of Model I 

In this appendix we describe in detail our first explicit MSSM orbifold compactification. 
We specify the orbifold, and present its complete 4D massless spectrum, along with the 
trilinear couplings allowed by the standard selection rules. Then we describe the modular 
symmetries and anomalies, and compute the coefficients that measure the stringy threshold 
corrections to the gauge kinetic functions. 

A.l Model definition 

The model is defined by the following shift vector and Wilson lines 

y so(io),i = fi _i _i )fi _I _i _i _i _I _i i) (A i) 

v y 3 , 2 ' 2' ' ' ' ' / V2' 6' 2> 2' 2' 2' 2 ' 2 I ' V^-^) 

M = (0, i, \, i, -i, 0, 0, 0) (1, 0, 0, 0, 2, -2, 0, 0) , 

A — (a — 1 —I —I —I —I -I 2 \ (10 _o _a _1 _1 _A _Q q\ 

\y, i, 3 , 3 , 3 , 3 , 3 , 3 j y 3 , <±, 3 , 3 , i, o, o j . 
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The Wilson lines .A3 (order 3) and A§ (order 2) lift the degeneracy of the fixed points in 
the second and third torus, respectively. A possible second order 2 Wilson line in the third 
torus is set to zero. The 4D gauge group after compactification is 

SU(3)c x SU(2) L x U(l)y x [SO(8) x SU(3) x U(l) 6 ] . (A.2) 

The U(l) generators are chosen to be 

ti = tx = (3, 1, —3, 3, 3, 3, 3, 3) (— 3, —3, 0, 3, 3, 0, 0, 0) , (A.3a) 

t 2 = t Y = (0, 0, 0, i, 1, -1,-1,-1) (0, 0, 0, 0, 0, 0, 0, ) , (A.3b) 

t 3 = (-45, 321, 45, -9, -9, -9, -9, -9) (9, 63, 0, -9, -9, 0, 0, 0) , (A.3c) 
U = (246, 0, 75, -15, -15, -15, -15, -15) (15, 105, 0, -15, -15, 0, 0, 0) ,(A.3d) 

t 6 = (0, 0, 171, 15, 15, 15, 15, 15) (-15, -105, 0, 15, 15, 0, 0, 0) , (A.3e) 

U = (0, 0, 0, 3, 3, 3, 3, 3) (54, -21, 0, -54, -54, 0, 0, 0) , (A.3f) 

t 7 = (0, 0, 0, 21, 21, 21, 21, 21 ) (0, 15, 0, 0, 0, 0, 0, ) , (A.3g) 

where ty and tx denote respectively the hypercharge and the anomalous U(l) generator. 

A.2 4D spectrum 

The spectrum in the untwisted sector includes: the graviton, the axio-dilaton, 5, the 
gauge bosons of the gauge group ( |A.2| ), the Kahler moduli Ti, i = 1,2,3, the complex 
structure modulus U3 (cf. Section fj~2]) , and their corresponding supersymmetric partners. 
In addition, the spectrum contains the matter fields with non-zero quantum numbers given 
in Table ||. They amount to three generations of quarks and leptons plus vector-like exotics 
with respect to Qsm and hidden matter. 

Table 3: The spectrum of the 2^-11 Model I. We denote the fixed points, /, in the 
•d k (un)twisted sector, with their (real) lattice coordinates. A pair of coordinates 
(— ,— ) indicates that the corresponding torus is left fixed in the sector. We give the 
representations under SU(3)c x SU(2)l x 50(8) x SU(3) respectively, and the U(l) 
charges qi with respect to the basis of generators given above. n l and £ l are the Kahler 
and complex structure modular weights in the ith torus. 
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A. 3 Trilinear couplings 

The string selection rules allow us to identify all non-vanishing couplings present in the 
effective theory emerging from an orbifold. Here, in Table ^j, we give all trilinear couplings 
allowed by the standard selection rules; gauge invariance, 7£-charge conservation and the 
space group selection rule. 
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Table 4: All couplings at trilinear level allowed by the standard selection rules. We display 
all 116 couplings, classified according to the twisted sector in which the involved fields 
appear. 



A. 4 Modular symmetries 

The target space modular symmetry for the first plane is the full group S'L(2,Z). The 
Wilson lines in the second and third planes break the symmetry in those directions down 
to congruence subgroups, which can be computed by solving the constraints Eqs. (48-55) 
found in the first reference of [25]. In the second plane the result is Ti(3). In the third 
plane, we have To(4) acting on the Kahler modulus, and r°(4) acting on the complex 
structure modulus. Here we collect some of the properties of these groups. 
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The fundamental domain for SL(2,Z) can be found in any text book on modular 
forms. It has a single cusp at infinity. Representing the modular transformations with the 
matrices 

( ab^ 
[ad 

they are generated by 



(A.4) 



. T: ^ Q1 j. (A.5) 

SX(2,Z) has two self-dual points, also known as elliptic points. They can be computed 
by solving = for some a,b,c,d and 0*. It is convenient to rotate <p so that its 

fundamental domain lies in the upper-half complex plane, taking — > —i<f>, in which case 
the elliptic point is given by: 





2c 

after using ad — be = 1. Here, c / and in order to lie within the upper half plane, we 
must have \a + d\ < 2. The solutions are at the boundary of the fundamental domain, 
at (returning to our original convention) = 1, ^ — |, which are fixed under S,TS 
respectively. 

Formally, the lowest weight modular form for SL(2, Z) has weight k = 12, and is called 
the discriminant function, A. It is related to the Dedekind eta function by A = r]((f)) 24: . 
There are also modular forms with so-called character, which transform up to a phase 
factor, and the Dedekind eta function with weight k = ^ is one of these. 

Several results for the congruence subgroups can be found in the textbooks. Since these 
groups are contained within S'L(2,Z), the above modular forms are also modular forms of 
To(N),T (N) and Ti(N). However, we also have additional ones, including r](N(/)) for 
Tq(N),Ti(N). Naturally, for a given N the largest space of modular forms will be for 
Ti(N). Generators for the groups can most easily be computed by using the computer 
program SAGE, although the resulting set is not small, nor even necessarily minimal. For 
Ti(3) we obtain a set of 12, of which we write down some examples: 

(A.7) 

It turns out that most congruence subgroups do not have elliptic points. Following 
the computation above, we find that this is so for To (4) and r°(4). For Ti(3) we find that 

</>* = ^ + § is fixed under ^ ^ ^ ^ . 
A.5 Modular anomalies 

The modular anomaly coefficients of this model can be computed with the modular weights 
n l , £ 3 provided in Table |3] by using Eq. ( 3.28| ). The result is given in Table |. Since all 
elements of the point group act non-trivially in the first torus, the anomaly coefficients b'^ 
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are universal. In contrast, they are different in the second and third planes for each of the 
groups. 
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Table 5: Non-Abelian modular anomaly coefficients b'i of a concrete Zq-II model. For 
j = 3, both Kahler-non-Abelian and the complex structure-non-Abelian mixed anomalies 
coincide, as expected [19]. 



A. 6 Auxiliary Af = 2 theories 

In Z^-II orbifolds the twist is given by v = |(0, 1, 2, —3). Taking the twist v% 2 = 3v 
and i> Z3 = 2v leads to two independent M = 2 theories. These are important for the 
computation of the stringy threshold corrections. Let us study each case separately. 

%2 theory. In the former case, v% 2 gives rise to a TL% orbifold with M = 2 in which the 
second torus is left invariant. Embedding the compactification into the gauge degrees of 
freedom amounts to taking the shift vector V% 2 = S^ 80 *- 10 )' 1 and the same Wilson lines 
A5, A3 as in the M = 1 orbifold (see Eq. ( AT)). Notice that the Wilson line A3 associated 
with the invariant torus acts non-trivially on physical states. 
The resulting M = 2 theory has the gauge group 

SU(5) x SO(10) x SU(3) x U(l) 5 . (A.8) 

Omitting the U(l) charges, the quantum numbers of the matter spectrum are given by 24 

18(5, 1, 1) 8(1, 1, 3) 4(1, 1, 1) c.c. (A.9) 

It is straightforward now to compute the beta function coefficients, which will be used 
below. The general formula is 

{V a f= 2 = -2C(G a )+Y J T{K), (ATO) 

a 

where i?^ denote the representations w.r.t. the non-abelian group G a , and the summation 
runs over all half-hypermutiplets of the theory. In the present case, for the non-Abelian 
groups we obtain 

(bl m )" =2 = 8, (6§o ( io)^ = -16. (^u(3)) A ' =2 = 2. (ATI) 



4 We list half-hypermutiplets. 
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7L<$ theory, yields a Z3 orbifold which preserves N = 2 and leaves the third torus 
untouched. The embedding into the gauge degrees of freedom requires 
and the same Wilson lines as before. 

The resulting N = 2 theory has the gauge group 

SU(3) a x SU(3) 6 x SO(8) x SU(3) C x U(l) 6 . (A.12) 

Omitting the U(l) charges, the matter spectrum is given by 

6 [(3, 1, 1, 1) © (1, 3, 1, 1) © (1, 1, 1, 3) © (1, 1, 1, 3)] © (3, 3, 1, 1) (A.13) 
© 3 [(3, 1, 1, 1) © (1, 3, 1, 1)] © 7(1, 1, 8, 1) © 17(1, 1, 1, 1) © c.c. 

Hence, the corresponding non-Abelian J\f = 2 beta function coefficients are given by 

(& 3 su(3) a f =2 = (bh(3 )b f =2 = (bh(3)J M=2 = 6 , (blo i8) f =2 = 2 • (A.14) 
A. 7 Universal S GS 



At this point, we have all ingredients to compute 5 GS as prescribed by Eq. ( 3.27] ) . Since 



our compactification has universal Kac-Moody level k a = 1, we must simply use 



For the second torus we have |P2|/|P| = 1/3 whereas for the third torus |P3|/|P| = 1/2. 
Furthermore, one has to be careful with how the 4D gauge group is embedded in the larger 
M = 2 theories. We find e.g. that SU(2)i is a subgroup of SU(5) in the %2 subsector, 
implying that for SU(2)^ 5 GS = &g 2 u( 2 ) i — l(^su(5))"^ =2 ' ^ ne can rea dily verify then that 

5gs = -Ti $GS = -h S 3 GS = 0. (A.16) 

Note that the same 5 GS enters the threshold correction contributions containing T3 and 
Us- Besides, 8 GS = holds always for Z^-II orbifolds 25 . 

B. Details of Model II 

In this appendix we describe in detail our second explicit example of an MSSM orbifold 
compactification, in an analogous way to that above. 

B.l Model definition 

The model is defined by the following shift vector and Wilson lines 

t/SO(10),1 _ (1 _1 _1 n n n n nWI —I —I —I —I —I —I l\ /r> i\ 

v 2 , 2 , u, u, u, u, uj y 2 i 6' 2' 2' 2' 2' 2' 2) ' \ D - L ) 

A = (—3 —I _I _1 _! 1 13 13W5 _3 _3 _3 _3 _3 _3 3 \ 
•S±5 \ 4) 4) 4) 4; 4 > 4 j 4 j 4 ) \ f ) 2' 2' 2' 2' 2' 2 ' 2 / ' 

^3 = {-\, -|, i i i |, f , f ) (0, 0, 0, 0, 0, 0, 0, 0) . 



D This is in general true in torii where the condition |P|/|Pj| = 2 is satisfied [19]. 
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The Wilson lines .A3 (order 3) and A§ (order 2) lift the degeneracy of the fixed points in 
the second and third torus, respectively. A possible second order 2 Wilson line in the third 
torus is set to zero. The 4D gauge group after compactification is 

SU(3) C x SU(2) L x U(l)y x [SO(14) x U(l) 5 ] . (B.2) 

The U(l) generators are chosen to be 

h = t x = (0, 0, |, HI |, |) (0, -2, 0, 0, 0, 0, 0, ) , (B.3a) 

t 2 =t Y = (0, 0, 0, -1,-1,1,1,1) (0, 0, 0, 0, 0, 0, 0, ) , (B.3b) 

t 3 = (3, 0, 0, 0, 0, 0, 0, 0) (0, 0, 0, 0, 0, 0, 0, 0) , (B.3c) 

t A = (0, 3, 0, 0, 0, 0, 0, 0) (0, 0, 0, 0, 0, 0, 0, 0) , (B.3d) 

t 5 = (0, 0, 42, -3, -3, -3, -3, -3) (0, 9, 0, 0, 0, 0, 0, 0) , (B.3e) 

te = (0, 0, 0, 9, 9, 9, 9, 9) (0, 15, 0, 0, 0, 0, 0, 0) . (B.3f) 

Again, ty and tx denote respectively the hypercharge and the anomalous U(l) generator. 

B.2 4D spectrum 

The spectrum in the untwisted sector includes: the graviton, the axio-dilaton, 5, the 
gauge bosons of the gauge group ( |B.2[ ), the Kahler moduli Tj, i = 1,2,3, the complex 
structure modulus Us (cf. Section |4~2| ), and their corresponding supersymmetric partners. 
In addition, the spectrum contains the matter fields with non-zero quantum numbers given 
in Table ||. They amount to three generations of quarks and leptons plus vector-like exotics 
with respect to Gsm and hidden matter. 

Table 6: The spectrum of the 1*6-11 Model II. We denote the fixed points, /, in the 
"9 k (un)twisted sector, with their (real) lattice coordinates. A pair of coordinates 
(— ,— ) indicates that the corresponding torus is left fixed in the sector. We give the 
representations under SU(3)c x SU(2)l x 50(14) respectively, and the U(l) charges 
qi with respect to the basis of generators given above. n l and t % are the Kahler and 
complex structure modular weights in the ith torus. 
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B.3 Trilinear couplings 

The string selection rules allow us to identify all non-vanishing couplings present in the 
resulting effective field theory. In Table ^ we show all couplings allowed by the standard 
selection rules; gauge invariance, 7^-charge conservation and the space group selection rule. 

B.4 Modular symmetries 



Much of the discussion is as in Appendix A. 4, The target space modular symmetry for the 



first plane is the full group SL(2, Z), whilst again, Wilson lines break the symmetry in the 
second and third planes. The groups can be computed to be Ti(3) for the second plane, 
Ti(2) acting on the Kahler modulus of the third plane, and T 1 (2) acting on its complex 



structure modulus. The fixed points of the first and second plane are as in [A.4j . For r^(2), 
transformation I I , leaves the point <fi* = 1 + i fixed. 



there is a fixed point at <^>* = \ — |, under the transformation ( ] . For T 1 (2), the 




B.5 Modular anomalies 

The modular anomaly coefficients of this model are computed with the modular weights 
n*, I? provided in Table [6| via Eq. ( |3.28 ). The result is given in Table ||. 



B.6 Auxiliary AT = 2 theories 

As in Model I, there are two N = 2 theories, which we study separately in the following. 
Z2 theory. The resulting N = 2 Z2 theory has the gauge group 

SU(5) x SO(16) x U(l) 4 . (B.4) 
Omitting the U(l) charges, the quantum numbers of the matter spectrum are given by 26 

18(5, 1)0 4(1, 1)0 c.c. (B.5) 

26 We list half-hypermutiplets. 
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Table 7: All trilinear level couplings allowed by the standard selection rules in Model II. 
We display all 160 couplings, classified according to the twisted sector in which the involved 
fields appear. 



j 

a 


1 


2 


3 


SU(3) C 


-11 


-1 


3 


SU(2) L 


-11 


-1 


3 


SO(14) 


-11 


-13 


-5 



Table 8: Non-Abelian modular anomaly coefficients b£ of Model II. For j = 3, both 
Kahler-non-Abelian and the complex structure-non-Abelian mixed anomalies coincide, as 
expected [19]. 
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The beta function coefficients are given in this case by 

(blv(B)f =2 = 8, (&io (16 )) Ar=2 = -28. (B.6) 
Z3 theory. The resulting ftf = 2 Z 3 theory has the gauge group 

SU(3) a x SU(3) b x SO(14) x U(l) 5 . (B.7) 
Omitting the U(l) charges, the matter spectrum is given by 

6[(3, 1,1)0(1, 3,1)] 0(3, 3,1) (B.8) 

e 3 [(1, 3, 1) e (3, 1, 1)] e 7(1, 1, 14) © 17(1, 1, 1) © c.c. 

Hence, the corresponding non-Abelian M = 2 beta function coefficients are given by 

(& 3 su(3) f =2 = (blv(3 )b f =2 = 6 , (bloiu))^ 2 = -10 • (B.9) 

B. 7 Universal S GS 

As for Model I, using the previous ingredients, it is straightforward to compute now S J GS . 
The result is 

6 GS = -11, 6 GS = -%, S 3 GS = 0. (B.10) 
Note that the same 5 GS enters the threshold correction contributions containing T3 and 

u 3 . 

C. Z 6 — II Trilinear (Twisted) Coupling Strengths 

In this appendix we compute the trilinear coupling strengths for twisted fields, for a generic 
^6-11 orbifold. These results can then be straightforwardly applied to both our explicit 
models. Then we describe how to show that the couplings are modular invariant. The 
allowed Yukawa couplings between twisted fields are of the kind t? 3 i9 4 ?9 5 and -# 2 i9 5 ?9 5 . Al- 
though we use only the latter in the main text, here we present them both for completeness. 

C.l Couplings #tf 4 tf 5 

Since i9 3 and i9 4 leave respectively the second and third torii fixed, J = {2, 3} in Eq. ( 3,12| ). 



The coupling strength is then given by 

y$ 7 = F a(3j9s VhN U5 ^exp{-2vrr 1 n t L 1 u}, (C.l) 

u 

where 

Ll= (-3 2 3 ) ' NU5 = 33/^3(1/3) (C - 2) 

and the 2D vector u reduces to the first two coordinates of u = fp — f a + A with A G A. We 
have computed all possible coupling strengths, and we find that there are only four different 
types of couplings. These are defined by the fixed points and thus are characterised by the 
difference f a p = fp — f a . We provide below the explicit forms and leading contributions. 
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1) /« = (0,0), 7/3 = (0,0), / 7 = (o,o). 

fop = (0,0). F aPl = 1 

Dominant contribution: exp{—2Tru t LiuTi} ~ 1 + . . . 
Zg~II sample coupling from Model I: 5\ n^Q Sj 

2) / Q = (i,0),/ /3 = (0,0), / 7 = (0,0). 
fap = (-1/2,0). F a ^ = V3 

Dominant contribution: exp{— 2nu t L\uTi} ~ e _7rT i/ 2 + . . . 
^6-11 sample coupling from Model I: 7130 5j 

3) / Q = (0,0), //3 = (i,0), / 7 = (0,0). 

= (1/3,0). F Q/37 = ^2 
Dominant contribution: exp{— 2HU 1 L\uTi} ~ e _7rT i/ 3 + . . . 
Zg-II sample coupling from Model I: <5i 7142 £7 

4) / a = (±,0), ^ = (±,0), / 7 = (0,0). 
faf) = (-1/6,0). F a/37 = ^6 

Dominant contribution: exp{— 27rti*Li?jTi} ~ e -71 " 7 !/ 3 + . . . 
Zg~II sample coupling from Model I: ^3 n$ 2 ^7 

C.2 Couplings tf 2 ?? 5 ?? 5 

Since t? 2 leaves the third torus invariant, J = {3} in Eq. ( 3.12| ). The coupling strength is 
then given by 

y Q 2 | 7 = F aM g s Vhh N 255 exp{-7r(Ti +T 2 u\L 2 u 2 )} (C.3) 

Jil,Ji2 



where Li is given in Eq. ( |C.2| ) and 

(2 -l\ 2 14 / 3 vr 4 

La = ^-1 2 J ' ^ = 3 3/2r6(i/3) (C - 4) 

and the 4D vector u = (ui,u 2 ) corresponds to the coordinates on the first two torii of 
u = fp — f a + A, where A is a lattice vector. There are again four types of independent 
couplings as given below, with their leading contributions. 

1) fa = fp = h = (0,0,/ su(3) ) with / su(3) = (0,0) or (§, §) or (§,±). 
f a p = (0,0,0,0). F aPl = 1 

Dominant contribution: exp{— vr(Ti u\L\U\ + T 2 u 2 L 2 u 2 )} ~ 1 + . . . 
Zg-II sample coupling from Model I: d\ dg u 2 

2) f a = (i,0,/ su(3) ), U = h = (0,0,/su ( 3)) with / su(3) = (0,0) or (§,§) or (f , ±). 
fafl = (-1/3,0,0,0). F a ^ = V2 

Dominant contribution: exp{— tt(Ti v^\L\U\ + T 2 u^Liu^)] m e - 27rT i/ 3 + . . . 
Zg-II sample coupling from Model I: d 2 d§ u 2 

3) f a = (0,0, ^, |), //3 = (0,0,0,0), / 7 = (0,0, |, |) or permutations thereof. 
f a p = (0,0,-1/3,-2/3). F aM = l 
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Dominant contribution: exp{— tt(T\ u\LiU\ + T2 U2L2U2)} ~ e 27rT2 / 3 + . . . 
^6-11 sample coupling from Model I: ng 7152 n-64 

4) /a = (|,0, i, |), ^ = (0,0,0,0), / 7 =(0,0,§,±) or 
U = (|,o, |, i), /£ = (0,0,0,0), / 7 = (0,0,i,|) 
/a/3 = (-1/3, 0,-1/3, -2/3), (-1/3, 0,-2/3, -1/3). F a ^ = 
Dominant contribution: exp{— 7r(Ti u\LiU\ + T2 it^-^^)} ~ e~ 27r ( Tl+T2 )/ 3 + . . . 
Zq-II sample coupling from Model I: 7123 n§2 ^64 

C.3 Modular covariance of Yukawas 

We now illustrate how to check that the Yukawa couplings transform covariantly under 
the modular transformations. Let us first focus on the i9 2 i? 5 $ 5 singlet couplings that we 
used in the main text for Model I. These are the couplings n^n^n^, n^zn^'n^ n^n^n^, 
^23^59^66) and they are of kind 3) and 4) above. We label as /13 and /14 their respective 
contributions to the holomorphic superpotential. Recalling that fields of the same charges 
and modular weights can transform amongst themselves under the modular symmetry, it 
turns out that previous couplings also transform amongst themselves. Let us now show 
this explicitly. 

Consider the SX(2,Z) transformations on T\. The required transformation property 
for the couplings is of the form (see Eq. ( |3.18| ) ) 



(icTi +d)M * (C.5) 





for some field- independent matrix M. Here, we have used that the sum over modular 
weights in the first plane is + + = — 2 for all the couplings. It is enough to show 
that the couplings transform well under the generators S,T of SL(2,Z). 

Take for example the S transformation on the first torus of a coupling of the third 

kind 

fc 3 (1/21,22) = ^255 Yl exp{-7rai^ui}exp{-7rT 2 t4L2U2)} (C.6) 

«1,U2 

where u\ = A. Concentrating on the T\ dependent exponential only, we now apply the 
Poisson resummation formula 



£ exp {-vr(A + ffP{\ + f) + 2xie\\ + /)} 
AeA 

= j V exp{-7r(A* + e)*p- 1 (A* + e)-2vri/*A*) (C.7) 

V A Vdet P x f^ A * L J 

to the sum over u\. Next split the sum over integers (mi, ni) G A* into a sum over integers 
(pi, ni) with mi = —3p± — 1, — 3pi, — 3p\ + 1. In this way, /^(l/Ti, T2) can be written as a 
linear combination of h^{T\,T2) and h>4 : (Ti,T2) above, with coefficients proportional to iT\ 
as required. 

The computation for the congruence subgroup Ti(3), which acts on T2 is similar, but 
note that one needs to apply the Poisson resummation formula twice. 
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We can make an analogous analysis for the Model II couplings. There we considered 
the couplings n-jn^n^, nnn§%n§$, n-jn^in^-, «i7^6i n 62, which are of kind 1) and 2) above. 
It turns out that these couplings transform into each other, with the computation following 
exactly as above. 
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